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A detailed description is given of an 
electronic stochastic analyzer for use 
with direct "real-time" measurements of 
the conditional distributions needed for 
a complete stochastic characterization of 
pulsating phenomena that can be repre- 
sented as random point processes. The 
measurement system described here is 
designed to reveal and quantify effects 
of pulse-to-pulse or phase-to-phase 
memory propagation. The unraveling of 
memory effects is required so that the 
physical basis for observed statistical 
properties of pulsating phenomena can 
be understood. The individual unique 
circuit components that comprise the 
system and the combinations of these 
components for various measurements, 
are thoroughly documented. The system 
has been applied to the measurement of 
pulsating partial discharges generated by 
applying alternating or constant voltage 
to a discharge gap. Examples are shown 
of data obtained for conditional and 
unconditional amplitude, time interval. 



and phase-of-occurrence distributions of 
partial-discharge pulses. The results 
unequivocally show the existence of 
significant memory effects as indicated, 
for example, by the observations that the 
most probable amplitudes and phases- 
of-occurrence of discharge pulses 
depend on the amplitudes and/or phases 
of the preceding pulses. Sources of error 
and fundamental limitations of the 
present measurement approach are 
analyzed. Possible extensions of the 
method arc also discussed. 
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1. Introduction 



There are many types of naturally occurring 
pulsating phenomena that have statistical proper- 
ties which have not yet been adequately explained. 
Included in this category of phenomena are certain 
types of nerve impulses, pulsating fluid flow and 
droplet formation, bursts of electromagnetic radia- 
tion from extraterrestrial sources, geological distur- 
bances such as earth tremors, and pulsating 
electrical discharges specifically considered in this 
work. These phenomena may exhibit complex 
chaotic behavior manifested by an apparent high 
degree of randomness in the time of occurrence 
and magnitude of the impulse events. For some 



phenomena, the complexity of the impulse behav- 
ior may, in part, be a consequence of memory 
propagation between successive events. In develop- 
ing a better understanding of the physical bases for 
pulsating phenomena, it is essential to assess the 
effects of memory propagation. 

In the case of pulsating partial-discharge phe- 
nomena, it has already been shown that effects of 
memory propagation are significant [1-4]. Partial- 
discharge (PD) phenomena are of special interest 
because they are types of localized electrical dis- 
charges that occur at defect sites in electrical insu- 
lation. Partial discharges are often the precursors 
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to insulation failure and represent undesirable 
electrical noise sources under some conditions. 
The detection of PD pulses has been used to assess 
insulation performance and integrity [5]. It is also 
known [6, 7] that PD phenomena exhibit stochastic 
properties that depend on the nature of the defect 
site such as characterized by the types of materials 
present as well as their geometrical configuration. 
Partial discharges also produce physical or chemi- 
cal changes in the characteristics of the defect sites 
(PD - induced aging) that in turn produce changes 
in the stochastic behavior of the discharge [8-11]. 

Efforts have been underway in numerous labora- 
tories to quantify statistically PD patterns using 
computer-assisted measurement and analysis 
techniques [12-22]. The incentive for this work has 
been the development of so called "smart" PD 
detectors that employ pattern recognition to help 
identify the type of defect at which the PD occurs, 
e.g., to distinguish between a cavity in solid insula- 
tion and a metal particle in liquid or gaseous insu- 
lation. Unfortunately, progress in the development 
of reliable automated methods for PD pattern 
recognition has been hampered by a failure to 
understand the physical mechanisms that deter- 
mine the stochastic properties of PD phenomena. 
In general, present computer assisted PD-measure- 
ment systems simply do not provide enough refined 
information about the stochastic properties of PD 
pulses for a meaningful analysis. 

The purpose of the present work is to describe a 
real-time stochastic analyzer that can be used to 
quantify the stochastic behavior of a train of 
electrical pulses that may or may not be correlated 
with a periodic time varying excitation source, e.g., 
a sinusoidal voltage. The instrument described 
here is an extended version of one that was used to 
investigate the stochastic behavior of pulsating 
negative-corona discharges generated by applying a 
constant voltage to a point-plane electrode gap 
[1, 24, 25]. In addition to the conditional pulse- 
amplitude and time-separation distributions that 
could be measured with the previous system, the 
present system also allows measurement of a set of 
phase-restricted pulse-amplitude and phase-of- 
occurrence distributions. This latter capability 
makes the instrument suitable for investigating the 
stochastic behavior of partial discharges generated 
using alternating voltages. The data acquired from 
this system provide immediate determinations of 
the existence of pulse-to-pulse or phase-to-phase 
memory propagation effects. 

The measurement system described here can be 
thought of as a type of electronic filter that is 



inserted between the impulse source and a 
computer-driven multichannel analyzer (MCA) in 
which data on the desired conditional or uncondi- 
tional pulse distributions are accumulated. The 
unique features of the circuitry of this filter are 
documented here in enough detail to allow replica- 
tion. The present system design incorporates 
standard commercially available nuclear-instru- 
mentation components, where possible, such as 
time-to-amplitude converters and linear pulse 
amplifiers. Although the present system can be 
employed to investigate any type of pulsating 
phenomenon that can be converted to electrical 
signals, it was designed primarily for the measure- 
ment of relatively stationary PD-pulse phenomena 
generated by a constant or low-frequency alternat- 
ing voltage. The system may not be well suited for 
investigations of impulses that have repetition rates 
much greater or less than the PD-phenomena 
considered here; and it will not perform well for 
phenomena that exhibit highly nonstationary 
behavior, i.e., phenomena for which the stochastic 
properties change rapidly with time. 

The range of phenomena to which the present 
system can be applied is considered and the 
system's inherent limitations and sources of error 
are analyzed. Extensions of the technique and 
alternative approaches that rely primarily on analy- 
sis using computer software are discussed. Exam- 
ples are presented of results obtained for partial 
discharges generated in a point-to-solid dielectric 
electrode gap. 

2. Definitions 

In this section we introduce the parameters that 
define the types of stochastic processes which can 
be investigated with the electronic measurement 
system described here. We also define the various 
conditional and unconditional distributions that 
are measured with this system and indicate how the 
measured distributions can be used to gain insight 
into the physical bases for the process under inves- 
tigation. 

2.1 Random Point Processes 

The types of pulsating phenomena to be consid- 
ered here are those that can be represented by a 
marked random point process as defined by Snyder 
[26]. In order to represent the phenomenon as a 
point process, the pulses must occur at discrete 
times that can be readily defined. In the case of a 
periodic time-varying excitation, the events of 
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interest must occur at discrete phases. This re- 
quires that an occurrence time (or phase) can be 
meaningfully associated with a particular property 
of the pulse such as its amplitude. Difficulties can 
be encountered in satisfying the criterion for a 
point process if, for example, there is significant 
variability in the shapes of the pulses or if there is 
the possibility that successive pulses can overlap or 
otherwise become indistinguishable. Ideally there 
should be a reasonable uniformity in pulse shapes 
and the mean spacing between pulses should be 
much greater than the pulse widths. The types of 
pulsating partial-discharge phenomena to which 
the present measurement system have been 
applied generally satisfy the requirements for a 
point process. 

It is also assumed that the point process can be 
marked by some property of the pulse such as its 
amplitude, width, shape parameters, or area under 
the pulse. In order to consider the mark as a prop- 
erty of the pulse measured with the present system, 
it is necessary that the mark be converted to an 
electrical signal with a voltage that is proportional 
to the "size" of the mark. For reasons previously 
discussed [1, 27], the partial-discharge pulse ampli- 
tude has been selected here as an appropriate 
mark which is a measure of discharge magnitude. 
Since PD can be detected by different methods, 
e.g., optical, acoustical, and electrical [5], it is 
necessary to convert the observed response to an 
electrical signal as is normally done for purposes of 
recording data. 

If the occurrence of the pulsating phenomenon 
is correlated with an externally controlled time- 
varying excitation process such as a chopped light 
beam or, as is sometimes the case for PD pulses, a 
sinusoidal voltage, then it may be more convenient 
to specify the time of pulse occurrence relative to 
the times of the external excitation processes. For 
PD-pulses generated with an alternating sinusoidal 
voltage, it is desirable to consider the phase-of- 
occurrence of a pulse as defined by the phase of the 
corresponding applied voltage at the time of 
PD-initiation. 

The point processes under consideration here 
are assumed to be random in the sense that both 
the times-of-occurrence and the marks can exhibit 
statistical variability, e.g., it is not possible to pre- 
dict precisely when a given pulse will occur or what 
its amplitude will be. For processes excited by a 
well-defined controllable periodic source, it is also 
possible to define point processes that are fixed in 
time or phase but exhibit statistical variability in 
the mark. Such a process might be, for example. 



the sum of the areas under all PD pulses that occur 
in a specified phase interval of the applied voltage. 
The sum could be recorded at a fixed phase imme- 
diately following the time lapse of the phase inter- 
val. The measurement system described in this 
work allows determination of such phase-restricted 
sums of pulse areas or amplitudes. 

2.2 Measurable Quantities for a PD Process 
(Random Variables) 

The type of marked random point process under 
consideration here is a stochastic process specified 
by a countable set of discrete random variables of 
which time-of-occurrence (or equivalently phase- 
of-occurrence) is one of the variables. In this 
section we define the sets of random variables that 
apply to the measurement of pulsating partial dis- 
charges generated either with a constant applied 
voltage (dc) or an alternating (sinusoidal) applied 
voltage (ac). 

2.2.1 Random Variables for a dc-Excited PD 

Process A diagrammatic representation of a de- 
generated PD process is shown in Fig. 1. As previ- 
ously discussed [1, 24, 25], this process can be 
specified by the finite set {qi,ti}„, i= 1,2, ..,, « 
where qi is the amplitude of the ith PD pulse (usu- 
ally expressed in units of picocoulombs) and ti is 
the time at which this pulse occurs. The measure- 
ment system described here records time separa- 
tions between successive PD events rather than 
actual occurrence times. It is therefore more con- 
venient to specify the process in terms of the set 
of random variables {^i, qi , A?,_i}„, i =2, ..., n , 
where /^ti-i=ti —ti-i is the time separation be- 
tween the (j-l)th and ith events. To satisfy the 
requirements for a point process, it is desirable 
that the mean duration of the PD events, as mea- 
sured, for example, by the pulse widths, 8ti, be 
much smaller than the mean time separation be- 
tween successive events, i.e., (A^,) l> (S^,) for all 
values of i. If all time intervals are recorded, the 
time-of-occurrence of any pulse can simply be 
determined from the sum 






(1) 



As will be seen from the discussion below, data on 
the time separations between events are needed to 
assess pulse-to-pulse memory propagation effects. 
If memory effects are important, the random 
variables associated with the amplitudes and time 
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separations of successive pulses are not indepen- 
dent, e.g., the amplitude of any given pulse can 
depend on the time separation between that pulse 
and the previous pulse. 
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Fig. 1. Diagrammatic representation of a marked random point 
process. Shown are pulses "marked" with amplitudes 
qj,i=n —2,n -1,..., that occur at discrete times /; with corre- 
sponding time separations At, . 



2.2.2 Random Variables for an ac-£xcited PD 
Process If the PD process is generated with an 
alternating voltage, it becomes more convenient to 
specify the phase-of-occurrence of the PD pulse 
rather than the time-of-occurrence. An example of 
an ac-generated PD process is shown by the 
diagram in Figs. 2a and 2b. The excitation voltage 
indicated in Fig. 2a is assumed to be sinusoidal and 
is given by 



Fa(0 = Kosin(wO = ^'osin[<^(0] , 



(2) 



where £u/27ris the frequency, (f>(t) = o}t the phase, 
and Vo is the amplitude. The individual PD events 
are specified by the set of random variables 






i = l,2,...,n', j =l,2,...,m , where 



qij" and qij' are the amplitudes of the i th pulses to 
appear respectively in theyth positive and negative 
half-cycles of the applied voltage and t^j/and 
^ij are their corresponding phases-of-occurrence. 
The phases <^j* are restricted by definition to lie 
within the interval (0, Itt) for arbitrary/ and are 
thus related the phase at time / by 
4>(t) = 4»f +27r(j-l). 

The amplitudes for qij' and qij are observed to be 
of opposite signs as indicated in Fig. 2b. In some 
cases, as previously explained [3, 4], the occurrence 
of positive and negative PD pulses may be phase 
shifted relative to the positive and negative half- 
cycles, e.g., it may be possible for negative pulses to 
occur before the zero-crossing where the voltage is 
still positive. This is phyiscally a consequence of a 
fluctuating phase lag in the local electric-field 
strength at the discharge site. The phase shift will 



be denoted here by 8(f>, and is arbitrarily adjusted 
to a value such that 

^if E( — 8^,Tr — 8<j>),^iJ e (tt — S<^, 27r — S<^) 

(3) 
for all values of i and j . 
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Fig. 2. Diagrammatic representation of an ac-excited partial 
discharge process: a) sinusoidal excitation voltage, b) phase- 
correlated PD pulses. 



As in the case of dc-generated PD pulses, it may 
also be useful to specify the phase differences 
between successive events within a half-cycle. 
These phase differences will be denoted by 
A4>i:i-i.i =<^/5 -^i~-tj wherein 2. 

For some types of ac-generated PD processes, 
especially those that occur in the presence of solid 
dielectric surfaces, it is valuable in assessing phase- 
to-phase memory propagation effects to know the 
accumulated PD charge associated with each half- 
cycle as defined by 



Qf = S qu 



(4) 



where the summation is over amplitudes of all 
pulses that occur within a given half-cycle as speci- 
fied by their phases-of-occurrence defined in 
Eq. (3). 

It is possible with the system described here, and 
sometimes necessary, to record the number of indi- 
vidual voltage cycles. This is necessary if an assess- 
ment is to be made of memory propagation that 
extends back beyond the previous cycle. For most 
types of measurements described here however. 
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this information is not recorded. If no attempt is 
made to record the number of a given cycle, then 
the subscript ; can be dropped from the specifica- 
tion of the random variables that define the 
stochastic process for ac-generated PD. In this 
case, the appropriate designation of the random 
variables is qf , ^f , A<f>ifi^i, and Q*. A failure to 
include the subscript on the variable Q ~ will imply 
by default that the sum given by Eq. (4) applies to 
the half-cycle immediately preceding that in which 
the variables q,^ , c^f , are measured. 

In performing the measurement of an ac-gener- 
ated PD process, it is assumed that the excitation 
voltage given by Eq. (2) is instrumentally filtered 
out or otherwise subtracted from the PD signals. 
This is necessary to ensure that the recorded pulse 
amplitude is a true measure of the discharge 
intensity. 

2.3 Measurable Conditional and Unconditional 
Distributions 

2.3.1 Unconditional Distributions Because the 
variables such as pulse amplitude and phase- 
of-occurrence that describe the pulsating 
phenomenon (PD process) of interest are random, 
they can only be specified quantitatively in terms of 
statistical probability distributions. The uncondi- 
tional probability distribution pi(x) (sometimes 
referred to as the probability density function [28]) 
for a random variable, f , is defined such thatp^Qc) 
dx is the probability that ^ will assume a value that 
lies in the interval ;c to j: -(- dx. Here ^ can be any of 
the random variables that were defined in the 
previous section. 

Consistent with our earlier work [1, 24], we shall 
adopt the abbreviated notation for distribution 
functions whereby p ((x) dx is replaced hypo(^) d^. 
Thus, for example, po(,qn) dq„ is the probability 
that the «th pulse has amplitude between q„ and 
q„ + dq„ . There is no ambiguity in this notation if it 
is understood that the symbol used to designate the 
value of a random variable is the same as that used 
to define the variable. 

The distribution po(^) is unconditional in the 
sense that it gives the probability that the random 
variable will have a particular value independent of 
the past history of the process, e.g., independent of 
values for random variables associated with previ- 
ous events. The random variables, as defined here, 
correspond to particular discrete events in time 
associated with the random point process, i.e., 
^ = qn, At„-i, <l>„ , where the subscript n is assigned 



to the n th event. In cases where the events are not 
actually counted by the measurement process, the 
distributions such as poiq,,) and po(At„) are 
assumed to apply to arbitrary n. If events are 
counted relative to a specified time, then n is 
assigned a value, e.g.,po(q2) is the amplitude dis- 
tribution of the second pulse to appear in an arbi- 
trary positive half-cycle of the excitation voltage. 

23.2 Conditional Distributions If memory 
effects are important in a pulsating process, then 
the probabilities that the random variables associ- 
ated with a particular event will have specific values 
depend in general on the values for these variables 
that were assumed by previous events. The proba- 
bility that the /th PD pulse will have values for 
amplitude and time separation that lie in the ranges 
q; to qj + dqj and A^j-i to Atj-i + d(Atj-i) can 
depend on values of all previous qi and Af,_i where 
i <j . The existence of memory effects can be 
established by the measurement of conditional 
probability distributions. The system to be de- 
scribed here allows measurement of a set of condi- 
tional distributions for such variables as pulse 
amplitude and phase-of-occurrence. 

The conditional distribution ;ji(^j|Afj-i) is de- 
fined such that pi(qj\Atj-i) is the probability that 
thej'th pulse has an amplitude in the range qj to 
qj + dqj if its time separation from the previous 
pulse has a fixed value Atj-i. With the system 
described here, it is also possible to measure higher 
order conditional distributions such as 
p2(qj\Atj-t,Atj-2), where p2(qi\Atj-i,Atj-2)dqj is 
the probability that the amplitude of thejth pulse is 
in the range qj to qj +dqj if both Atj-i and Atj-2 are 
fixed. Lists of the conditional distributions that can 
be measured for dc and ac-generated PD pulses are 
given respectively in Tables 1 and 2. 



Table 1. Measurable conditional and unconditional pulse-ampli- 
tude and time-separation distributions for a constant excitation 
process (de-generated PD) 



Distribution 
type 


Amplitude 
distribution 




Time-separation 
distribution 


Unconditional 


Pofe) 




Pa(AO) 


First-order 
conditional 


pi(qj\Mj-i) 

Pifekj-i) 
pi(qi\Atj-;,),k 


>1 


PiiMjlMj-i) 


Second-order 
conditional 


P2(qj\Mj-i,qi- 
p2(qj\Atj-i.Ati. 


i) 

-2) 
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Table 2. Measurable eonditional and uneonditional pulse-amplitude and phase distributions for a periodie, time vaiying excitation 
process (ac-generated PD) 



Distribution 
type 


Amplitude 
distribution 




Phase 
distribution 


Pulse number and 
cycle specifications 


Unconditional 
(unspeeified event) 


P<i(qf ) 




po(<^f) 


all, (SI 


Uneonditional 
(specified event) 


Paiqf ) 




Po(<^f) 


(■ = 1,2,3,... 


First-order 
eonditional 


Pi{qf\<i>f) 

pAqf, \Qt ) 




Pi(<^f|<^f-.) 
Pi(<^/;lGf) 


<i>f eA<^*,all,i2:l 

( = 1,2,3... 

( = 1,2,3... 

+ -^k. =],]-!,... 

-^k=j-l,j-2,... 

i =2,3,4,... 


Seeond-order 
conditional 


P2(qf\<t>f-i,qr- 

P2{qf,\<^f,,Qk) 


.) 


P2(<l>^\<l>^-uqr-^) 


i =2,3,4,... 
or all (■ 5 2 

(=1,2,3,... 

+ ^k =j,j-l,... 

-^k=i-\,j-2,... 


Third-order 
conditional 


p:i{qf \<i>f ,qT-v 


.<^f-i) 


p:i(<l>fi\<l>f-\.i<qf-\.i>Qk ) 


i =2,3,4... 

+ -*k =J,J-l,... 

-^k=j-l,j-Z... 



Determination of the conditional distributions such 
as pi(qj\Atj-i) provides an indication of the de- 
pendence of the random variable Qj on Atj-i. If 
memory effects are important, then the probability 
that qj will assume a particular value can depend 
on the value chosen for Atj-i. In this case, the con- 
ditional distribution, pi(qj\Atj-i), will not equal 
the unconditional distribution, /?o(9;), for at least 
some allowed values of Atj-i. 

A quantitative assessment of memory propaga- 
tion can be made from calculation of expectation 
values using measured conditional distributions. 
For example, the expectation value for the phase- 
of-occurrence of the third pulse in a negative half- 
cycle of the applied excitation voltage conditioned 
on a fixed value for the sum of all PD pulse ampli- 
tudes in the previous positive half-cycle is defined 
by 



{<t>nQ^))= <t>ipi(.<t>i\Q^)A<i>i 

J-ir-SA 



(5) 



where it is assumed that <f>2. must be confined to 
the interval defined by Eq. (3). In general, 



<^.(W}'.))= f ^iP^{^i\{ak}n)A^i, (6) 

JR 

where ^/ is any random variable associated with the 
/th pulse and {a*}„ is a set of fixed values for n 
random variables associated with one or more 
pulses that occurred at earlier times. The integral 
in Eq. (6) is over all allowed values of ^, that are 
assumed to lie within a range R , i.e., ^,- € R. 

If memory effects are important, the value of 
{Si ({a* }h )) will change as one or more of the values 
flk are changed. If the value of (^, ({a*},, )) increases 
as aie{ak}n increases within a particular range 
(at €.Ai), then ^, is said to hepositively dependent on 
a I in that range for fixed values of a* (k^^l). Consis- 
tent with our earlier notation [1], this dependence is 
denoted by (a/ 1 ^ ^/ ti {ok,k ^l}) when a; €Ai. 
Likewise ^j can be negatively dependent on a differ- 
ent variable or on at in a different range ^4 /. In this 
case, the negative dependence is denoted by 
(a/ 1 => ^1 i, {okjk Til}) when at eA'i. The depen- 
dence of the expectation values for random 
variables on the values of random variables associ- 
ated with prior events can often be predicted from 
physical models of the process as has been done for 
the case of negative corona (Trichel) type partial- 
discharge pulses [1]. 
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If memory effects are important in the pulsating 
phenomenon, then the various distributions listed 
in Tables 1 and 2 are not necessarily independent. 
It can be shown, for example, from the law of prob- 
abilities that the distributions po(9j),Pu(A(j-i), and 
pi(qj\A.tj-i) are related by the integral expression 



All conditional distributions satisfy the normal- 
ization requirement 



JR 



(10) 



Pofe) = f PB(Ar;-i)Pl(g;|Afy-i)d(Af;-,). (7) 

Jo 

Similarly the distributions po(<^/~ ), PoiQ ~), and 
Pi{4>t \Q *) are related by 



Measured data for conditional distributions are 
generally normalized according to this requirement 
by numerical integration. In some cases, it may be 
convenient for display purposes to normalize to the 
maximum of the distribution. 



M«^f ) = f 



Po(Q")pi(<^f|e')dQ = 



(8) 



It has previously been shown [1] that there may 
be many other integral expressions that connect 
the different conditional and unconditional 
measurable distributions. Equations (7) and (8) 
can be used to check the consistency among the 
various measured distributions. For 

example, if data are obtained on the three distribu- 
tions pofe),P()(A/'j-i), and pi{qi\Ltj-i), then one 
should, if possible, verify that they satisfy Eq. (7). 
There may be some cases, however, where it is not 
possible to obtain enough data at high enough 
resolution to perform this analysis. 

In the process of measuring a conditional distri- 
bution it is generally not possible to select a single 
value for the "fixed" variable. This variable can 
only be specified experimentally to lie within a 
finite window. In the case of the distribution 
Pi(4>r \Q*), for example, one really measures an 
approximation to this distribution given by 



Jo* -so* 



(9) 



where Q'* is defined by the measurement to lie 
within the window corresponding to the interval 
(Q^ -SQ^,Q^ + SQ^). The measured conditional 
distribution approaches the "true" conditional 
distribution as the window is made smaller, e.g.. 



Pi(<l>r\Q^)= lim pi(<l>r\Q*±5Q*). 

SQ*-*0 

The errors associated with finite window size have 
previously been noted [24] and will be discussed 
again later in this work. 



3. Measurement System 

In this section we describe the general features 
of the system for measuring the conditional and 
unconditional distributions listed in Tables 1 and 2. 
The system can be configured to investigate either 
a continuous train of pulses produced by a constant 
excitation process, e.g., dc-generated PD or pulses 
generated by periodic, time-varying process, e.g., 
ac-generated PD. Thus, the system is an extended 
version of that previously described for measure- 
ment of dc-generated PD [1, 24, 25], and, in fact, 
includes all of the features of the earlier system. 
We shall treat the ac and dc measurement configu- 
rations separately even though they both utilize 
some of the same individual circuit components. 

3.1 Conflgurations for a Continuous Excitation 
Process (de-Generated PD) 

The configurations of the electronic system used 
to measure the distributions listed in Table 1 for a 
dc-generated PD process have been described 
previously [1, 24, 25], a block diagram indicating 
the circuit components utilized for this case is 
shown in Fig. 3. The configurations of the compo- 
nents that are required for the measurement of 
each distribution are specified by the various 
switch configurations listed in Table 3 where 
51-S7 are the switches designated in Fig. 3 and the 
notation a 1=61 implies, for example, that position 
fl 1 of SI is connected to position M of 51. 

The system indicated in Fig. 3 differs from that 
used previously mainly in the design of the individ- 
ual circuit components that will be described in the 
next section. The most significant changes have 
been in the design of the At control logic circuits 
(parts A and B). Although their basic function and 
operation are the same as previously described 
[1, 24], changes were made that reduce errors, 
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Fig. 3. System for measuring uneonditional and conditional pulse-amplitude and time-separation distributions for a continuous or dc 
exeited point proeess sueh as shown in Fig. 1. A = amplifier, DDG = digital delay generator, TAC = time-to-amplitude converter, 
SCA = single-ehannel analyzer, MCA = multichannel analyzer, G = gate, S = switch. 



Table 3. Configuration of switch conneclions for the system shown in Fig. 3 that are required for measurement of the various condi- 
tional and unconditional pulse-amplitude and time-interval distributions for a constant excitation process 











Switch 








Distribution 


SI 


S2 


S3 


S4 


S5 


S6 


S7 


Piife) 


xl=zl 


* 


« 


x4=y4 




* 


x7=z7 


Pi{cij A'j-i) 


xl=zl 


x2=y2 


x3 = w3 


x4 = x4 




* 


x7=z7 


Pi{qj Mj-t),k>l 


xl=zl 


x2=y2 


x3 = w3 


x4 = z4 




* 


x7=z7 


piiqj Mj-uqj-O 


xl=zl 


x2 = z2 


x3 = z3 


x4=w4 




* 


x7=z7 


P2(.qj Atj-u^tj-2) 


xl=zl 


x2 = z2 


x3=y3 


x4 = w4 




* 


x7=z7 


Po(Mj) 


xl=yl 


* 


« 


* 




x6 = z6 


x7=y7 


p,{M/\Mi-,) 


xl=yl 


* 


« 


* 


x5 = z5 


x6=y6 


x7=y7 


pi(.mqj) 


xl=yl 


» 


* 


* 


x5=y5 


x6=y6 


x7=y7 



* Switch position irrelevant. 



improve performance, and eliminate redundancies. 
Some of these changes have already been utilized 
in the investigation of Trichel pulses [1], but have 
never been documented in detail. 

The pulse sorter circuit that drives the time-to- 
ampUtude converters (TACs) and the analog gate 
G3 that accepts the outputs of the TACs remain 



unchanged. The TACs, digital-delay generator 
(DDG) and the 256-channel multichannel analyzer 
(MCA) are commercially available instruments. 
The commercial single-channel analyzers (SCAs) 
have been replaced with a circuit to be described in 
the next section. 
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The operating principles of the system shown in 
Fig. 3 have been discussed already in previous pub- 
lications. This configuration allows the recording of 
either pulse amplitudes or pulse time separations 
with a computer controlled MCA. The MCA 
employs a fast analog-to-digital converter to digi- 
tize the voltage amplitude of pulses received at its 
input provided the amplitude is above a preset 
discrimination level. If the MCA is set by switch S7 
to measure pulse amplitudes, then the input signals 
are derived from a linear pulse amplifier, Al, after 
passing through the analog gate G2. This gate is a 
built-in feature of the MCA. If time intervals are 
recorded, then the MCA input is derived from the 
output of one or more TACs. The output of the 
TAC is a narrow pulse with an amplitude that is 
directly proportional to the time between the start 
and stop pulses that are applied respectively to the 
"start" and "stop" inputs. For the measurement of 
the unconditional amplitude and time separation 
distributions (po(^y) and po^Atj)), the gate G2 is 
held continuously open by positioning switch S4 
such that x4=y4. 

The measurement of conditional pulse-ampli- 
tude distributions requires the use of the Ar control 
logic circuits and the digital-delay generator 
DDG3. For the measurement of the first-order 
conditional distribution /Ji(^;|Ar;-i), the gate, G2, 
to the MCA is enabled by the output of the Ar 
control logic (part A) for a time interval 
A^y-i±5(A^;-i) after the;- 1 th pulse, provided no 
pulse from Al occurs within the interval 25(Ar/_i) 
starting from the;-! th pulse. The time delay, Ar/_i 
and window ±5(Ary_i) are determined respec- 
tively by the settings of the DDG3 delay and pulse 
width. The ;-l th pulse essentially triggers DDG3 
after passing through the Ar control logic (part A). 
The DDG3 then returns a 5 V logic pulse of preset 
delay and width which is in turn transferred to G2 
if no other pulses have appeared at the input /. 
Details of the A^ control logic circuit operation are 
given in the next section. 

The measurement of the distributions 
pi(qj\Atj-i), where i>l and A^;-, is the interval 
between the; -/ th and; —i + lth events, requires 
the use of both parts A and B of the Ar control 
logic circuit. The value of / is determined by a 
selectable pulse counter inherent to the At control 
logic (part B) as described in the next section. The 
time interval Atj-i is selected by part A of the At 
control logic in conjunction with DDG3 as in the 
case of the pi(qj\Atj-i) measurement. The output 
of part A is then used to trigger part B. If a pulse 
appears at the input eb of part B within the time 



interval Atj-i ±S(Atj-i), then the gate G2 to the 
MCA is enabled either immediately for measure- 
ment o{pi(qj\Atj-2) or after i —2 pulses have been 
counted for measurement of pi(qj\Atj-i), i>2. 
The next pulse to appear after G2 is enabled will 
be recorded by the MCA which then returns an 
"event pulse" to reset part B of the At control 
logic. The next input pulse time interval to lie 
within the range Ary-,±5(Ar/-,) will start the 
process over again. 

The second-order conditional pulse-amplitude dis- 
tributions, /J2(9y |Ar;_i,^y_i) and p2(qj |Ar;_i,Ary_2) 
can be measured by using a single-channel ana- 
lyzer, SCAl, the output of which is connected to 
input e" of the Ar control logic circuit (Part A). 
Depending on the position of switch SI, SCAl 
receives a pulse either directly from amplifier Al 
or from the output of TACl for measurement of 
P2iqj\Atj-i,qj-i) and p2iqj\Atj-i,Atj-2) respec- 
tively. In the first case, the A^ control logic and 
subsequently DDG3 are only triggered if the amp- 
litude of qj-i lies within a narrow range selec- 
ted by SCAl. In the second case, it is triggered only 
if the output pulse of TACl, the amplitude of 
which is proportional to Ary-2, lies within a narrow 
range corresponding to Atj-2±S(Atj-2) as selected 
again by SCAl. 

The first-order distribution ;ji(^y|^y-i) can be 
measured using the configuration for measurement 
of;j2(^y |Ary_i,^y_i) and selecting the time window 
d(Atj-i) from DDG3 to be large compared to the 
mean time separation between pulses, i.e., 
d(Atj-i) > (Atj-i). Although it is possible to mea- 
sure directly other types of conditional amplitude 
distributions with this system [24] such as 
/Ji(^y|A^y-i -t-A^y_2), these are derivable from the 
distributions listed in Table 1, hence are consid- 
ered difficult to interpret and less useful in reveal- 
ing stochastic properties of the process. 

Measurement of the unconditional time-separa- 
tion distribution, ;jo(Ary), requires use of two time- 
to-amplitude converters (TACl and TAC2) 
connected to a pulse sorter. As previously shown, 
[24] this arrangement allows measurement of all 
successive time separations if all the time separa- 
tions are greater than the TAC reset time. The re- 
set time for the TACs used in the present 
measurement system is 50 |xs. Failure to sample all 
time separations can lead to errors in the measure- 
ment o{po(Atj) under some conditions as will be 
discussed later [24]. 

Measurement of the conditional time-separation 
distributions involves use of the single-channel 
analyzer SCA2 that is either connected at S5 
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directly to Al for measurement ofpi(Atj\qj) or to 
the output of TACl for measurement of 
pi(Atj\Atj-i). The output of SCA2 enables gate Gl 
for measurement of A/j with TAC3 provided either 
qj or Ary-i lie within the windows selected by 
SCA2. The gate Gl is actually a built-in feature of 
the time-to-amplitude converter circuit used in the 
present system. 

3.2 Conflguration for a Periodic Time-Varying 
Excitation Process (ac-Generated PD)] 

A diagram of the system configuration used for 
measurement of the distributions listed in Table 2 
is shown in Fig. 4. Although it is assumed here for 
convenience and simplicity that the excitation 
process for the observed pulses is sinusoidal as 
given by Eq. (2), this is not a requirement for the 
measurement method. It is only required that the 
excitation process have a well defined periodicity 



so that phase position and intervals can be mean- 
ingfully specified. Thus, excitation processes that 
can be represented by a Fourier expansion are also 
acceptable, e.g., voltages of the form 



Fa(0 = Co+2^n cos(n(ot)+^B„ sin(ncot), 

n=l n=l 

(11) 

where Co is a constant and A„ and B„ are the usual 
Fourier expansion coefficients. 

One of the major differences between the 
measurement system shown in Fig. 3 and that 
shown in Fig. 4 is that, in the latter configuration, 
the measurement of pulse occurrence times are 
always made relative to a fixed reference phase. 
This reference is provided by the output of a zero- 
crossing detector similar to that used in our earlier 



AC 

input 

O H 



Zero 
crossing 
detector 



Pulse counter 
cycle selector 




b" 



Id" 



SCA1 



Absolute 
value (+or-) 



Gate 



S1 



t-^ 



DDG1 



DDG4 




Gated 
integrator 



f 



Pulse 
selector 



x1 



y1?z1 
+5V 



e" 



z5 

w5 



c"S2^DDG2 



x2 9y2 



b" 



SCA2 



b" 



Start 



Stop 



TACl 



*->■ 




A3 



6S5 



y5 



O x5 e' 



At 

control 

logic 



w3l 



DDG3 



f" 




+5V y6 , _ 

CMD 1 6z6 
366x6 



Fig. 4. System for measuring unconditional and conditional pulse-amplitude, amplitude-sum, and pliase-of-occurrence distributions for 
a point process excited by an alternating voltage. The individual circuit units are defined as in Fig. 3. 
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work [29] which generates a positive 5 V pulse with 
a width of 2 fjus whenever the excitation voltage 
changes sign from negative to positive. The output 
of the zero-crossing detector is fed to a pulse coun- 
ter. The output of the pulse counter triggers two 
digital-delay generators DDGl and DDG4 used to 
define the phase intervals over which measurements 
are made. Depending on the setting of the pulse- 
counter output, DDGl and DDG4 are triggered 
either at the beginning of every cycle of the excita- 
tion voltage or at the beginning of every nth cycle, 
where n is an integer greater than 1. 

The MCA, TAG, gate Gl, At control logic (part 
A), SCAs and DDG3 are the same components 
used for the system configuration shown in Fig. 3. 
The other digital-delay generators DDGl, DDG2, 
and DDG4 are essentially identical in their operat- 
ing characteristics to DDG3. The gated integrator 
and pulse selector are specifically designed for the 
system configuration shown in Fig. 4 and are 
described in the next section. The absolute-value 



selector circuit is similar in design to a circuit used 
for previous PD measurements [30]. It provides a 
positive pulse to amplifier A2 independent of the 
sign of the input pulse. It can be operated to select 
either positive input pulses, negative input pulses, 
or pulses of both signs. This feature is needed be- 
cause ac-generated PD pulses are either positive or 
negative depending on the half-cycle of the applied 
voltage in which they occur. The amplifier A2 is a 
commercial linear pulse amplifier that has a con- 
stant adjustable dc offset at the output. It delivers a 
rectangular negative pulse to the gated integrator 
with a constant width of 2 fjus and with an ampli- 
tude proportional to the peak amplitude of the 
input pulse. It has output characteristics required 
for proper operation of the integrator. Shown in 
Tables 4 and 5 are the combinations of switch 
connections for S1-S5 in Fig. 4 that are required to 
configure the system for measurement of the 
various conditional and unconditional distributions 
given in Table 2. 



Table 4. Configuration of switch connections for the system shown in Fig. 4 required for measurement of the various conditional and 
unconditional amplitude or total charge distributions for a periodic time-varying excitation process 









Switch 








Distribution 


SI 


S2 


S3 


S4 


S5 


S6 


po(?*),/>l 


* 


* 


x3 = y3 


x4 = z4 


« 


x6 = z6 


Pofef) 


xl = zl 


x2 = y2 


x3 = w3 


x4 = y4 


x5 = z5 




pa(qr),i=2,3,... 


xl = zl 


x2 = y2 


x3 = w3 


x4 = y4 


x5 = w5 




Pi(9r|A<^*) 


« 


* 


x3 = w3 


x4 = z4 


« 


x6 = z6 


P>iir\'t>r) 


xl = zl 


x2 = y2 


x3=w3 


x4=w4 


x5 = z5 




P,{qr\<i>r\i=^X... 


xl = zl 


x2 = y2 


x3 = w3 


x4 = w4 


x5 = w5 




Pi{qt\Q^) 


xl = zl 


x2 = z2 


x3 = w3 


x4 = y4 


x5 = z5 




Pi(9*|G '),'■ = 2,3,... 


xl = zl 


x2 = z2 


x3 = w3 


x4 = y4 


x5 = w5 




P2{qf\<i>r,Q') 


xl = zl 


x2 = z2 


x3 = w3 


x4 = y4 


x5 = z5 






xl = zl 


x2 = z2 


x3 = w3 


x4 = y4 


x5 = z5 




P2(qr\<t>r,Q^),i =2,3,... 


xl = zl 


x2 = z2 


x3 = w3 


x4 = w4 


x5=w5 




p,(9r|A<^i:,-,)-'=2,3,... 


xl=zl 


x2 = y2 


x3 = w3 


x4 = y4 


x5 = w5 




P2(9/*|<^,li.?r-i),(^2 


« 


* 


x3 = w3 


x4 = y4 


x5 = y5 




P3(?;-|<^*,<^f-l?/ll),/^2 


* 


* 


x3 = w3 


x4 = y4 


x5 = y5 




P»iQn 


« 


« 


x3 = z3 


x4 = z4 


• 


x6 = y6 



* Switch position irrelevant. 

Table 5. Configuration of switch connections for the system shown in Fig. 4 required for measurement of conditional phase distribu- 
tions for a periodic, time-varying excitation process 











Switch 








Distribution 


SI 


S2 


S3 




S4 


S5 


S6 


P„(<^r),' = i,2... 


xl = zl 


x2=y2 


x3=y3 




x4 = z4 


« 


x6 = y6 


Pi(<^f|G '),'■ = 1,2,... 


xl = zl 


x2 = z2 


x3 = y3 




x4 = z4 


* 


x6 = y6 


Pi(<t>r\<t>r-i.),i=2,3... 


xl=yl 


x2 = y2 


x3 = y3 




x4 = x4 


x5 = y5 


x6 = y6 


P2(<^f|<^f-i,9;-i),'=2,3... 


xl=yl 


x2 = y2 


x3 = y3 




x4 = z4 


x5 = y5 


x6 = y6 


Pi(<l>r\<t>r-i,qr-i,Q^),i = 2,3,.. 


. xl = yl 


x2 = z2 


x3=y3 




x4 = z4 


x5 = y5 


x6=y6 



* Switch position irrelevant. 
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An understanding of the system operation can 
be obtained with the aid of the pulse diagrams 
shown in Figs. 5-8. Consider first the measurement 
of the unconditional amplitude distribution, 
Po{qf), of the ith pulse in a particular half-cycle. 
Figure 5 shows the time sequence of signals that 
appear at the various indicated circuit locations in 
the system shown in Fig. 4 for the measurement of 
po{q \), i.e., the first pulse to appear on the nega- 
tive half-cycle. 



Applied 
voltage 



X2.C' 




Fig. 5. Fulse timing diagram for indicated signal locations in the 
system configuration (Fig. 4) for measurement of the amplitude 
distribution, /7o(9 r)> of ^^^ fi"^' pulse to appear on the negative 
half-cycle of the applied voltage. 
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Fig. 6. Pulse timing diagram for indicated signal locations in the 
system configuration (Fig. 4) for measurement of the condi- 
tional distribution pi((fc*|2"). 
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Fig. 7. Pulse timing diagram for indicated signal locations in the 
system configuration (Fig. 4) for measurement of the condi- 
tional distribution P2(9 i\<t>I<Q*)- 
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Fig. 8. Pulse timing diagram for indicated signal locations in the 
system configuration (Fig. 4) for measurement of the condi- 
tional distribution p2(<^ 2 I <^r>9 D- 

The cycle selector is assumed here to be set such 
that every pulse a ' from the zero-crossing detector 
triggers the digital-delay generator DDGl. The 
output pulse from DDGl is delayed relative to the 
zero-crossing pulse and its width is adjusted to 
encompass the entire phase region within which 
the negative pulses occur. This pulse is fed through 
switch S2 to the pulse-selector circuit which also 
receives pulses at location f from the pulse ampli- 
fier plus absolute value circuit. The pulse selector 
always generates a "start" pulse at e which is 
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coincident with the leading edge of the DDGl 
pulse at a fixed phase <f>s. The start pulse triggers 
the At control logic circuit (part A) at e' which in 
turn triggers DDG3. The position and width of the 
DDG3 output pulse is adjusted in this case to be 
approximately the same as the DDGl pulse. The 
output from the At control logic opens gate G2 to 
the MCA. It disables this gate after the next pulse 
appears at the f input (in this case the first negative 
pulse at <f> i). The amplitude of the first negative 
PD pulse is then recorded by the MCA shortly 
before G2 closes. 

Measurement of po(qf) for values of i greater 
than 1 requires that the At control logic be trig- 
gered by the (/ — l)th pulse. This is achieved by 
connecting e' to e through switch S5 (x5 =z5). The 
pulse selector will produce a pulse at e coincident 
in time (or equivalently in phase) with a 
selected pulse that occurs within the range defined 
by the DDGl pulse. For example, 'ifpo(q J) is to be 
measured, then the second pulse is selected to 
appear at e. This pulse triggers the At control logic 
which in turn allows passage of the third pulse to 
the MCA. 

The measurement of phase-of-occurrence distri- 
butions requires that the output of the time-to- 
amplitude converter, TACl, be recorded by the 
MCA. The TACl circuit is triggered by the outputs 
from the pulse-selector circuit. The diagram in 
Fig. 6 shows the sequence of signals at the indi- 
cated circuit locations associated with the measure- 
ment of the conditional phase distribution 
Pi(<f>2\Q~)- Again, it is assumed that the cycle 
selector allows every zero-crossing detector pulse 
to trigger DDGl. In this case, the output pulse 
from DDGl is used to control a gated integrator. 
The width and position of the DDGl pulse is set to 
encompass all possible negative pulses that could 
occur in the negative half-cycle. 

The integrator returns a pulse at the end of the 
DDGl pulse that has an amplitude proportional to 
the sum, "Zqr , of the amplitudes of all negative 
pulses contained within the phase window defined 
by the DDGl pulse. If this sum has a value that lies 
within the window Q~ ±8Q~ defined by the single- 
channel analyzer SCA2, then SCA2 triggers DDG2 
which controls the pulse selector. The pulse selec- 
tor produces a TAC start pulse at the leading edge 
of the DDG2 pulse at a fixed phase </>s. For the 
example shown in Fig. 6, the pulse selector is set to 
select the second positive pulse to appear in the 
window defined by DDG2. This pulse is used to 
stop TACl which then produces a pulse of ampli- 
tude directly proportional to A</>2 = <f>2—(f>s- The 



actual phase-of-occurrence of the second pulse 
relative to the zero crossing is given by 



c{>t-- 



■■({>s-2Tr + KArAc, 



(12) 



where Atac is the amplitude of the TAC output 
pulse and k is a scale factor determined from a 
calibration of the TAC. By this process, <f>2 is 
recorded only if the sum of negative pulse ampli- 
tudes in the previous half-cycle lie within a re- 
stricted range thus yielding the conditional 
distribution pi(<t>3\Q~). 

Figure 7 shows the sequence of signals associ- 
ated with the measurement of the second-order 
amplitude distribution, /72(9r|</> 7, 2^). As in the 
case of the po(q J") measurement indicated by the 
diagram in Fig. 5, this measurement requires use of 
both the pulse selector and At control logic circuits. 
Unlike the measurement of po(tj F), the amplitude, 
^ r, is recorded only if the phase of this pulse and 
the sum of amplitudes on the previous half-cycle 
have values that lie within restricted ranges. The 
phase, <f>T, is restricted to lie within a range speci- 
fied by the delay and width of the DDG3 pulse. As 
in the case for measurement of/7i(</>f |j3''),j3'' is 
specified by the SCA2 window. 

There are two switch configurations that can be 
used to measure p2{qf\(f>f ,Q*). The diagram in 
Fig. 7 corresponds to the first set of switch connec- 
tions for this distribution listed in Table 4. The 
second set of switch connections are required for 
measurement of P2(qf \<f>f ,Q*) if i>2. The 
phase, <f>f , in this case is restricted not by DDG3, 
but rather by DDG4 that controls gate G3. The gate 
G3 is actually a built-in feature of the At control 
logic circuit (part A) as will be shown later. 

The case for/72(</>2 l</>r,9r) is shown in Fig. 8. 
This measurement also requires use of both the 
pulse selector and At control logic circuits. The At 
control logic is used to enable a gate in the pulse 
selector that controls passage of the TAC stop 
pulse. The amplitude q f is restricted by the window 
setting for SCAl. The phase, </>f, is restricted by 
the DDG4 pulse which controls SCAl. If q f lies 
within the specified range, SCAl will produce a 
pulse at the end of the DDG4 pulse that triggers the 
At control logic. The At control logic uses DDG3 to 
open the pulse selector gate. If the next pulse to 
occur is indeed the second pulse, it will cause a stop 
pulse to be transmitted to TACl. If the next pulse 
is not the second pulse, then the stop pulse output 
of the pulse selector will simply be disabled. Thus 
only the phase of <f>2 can be recorded provided the 
previous pulse satisfies the specified conditions for 
amplitude and phase. 
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The measurement of unconditional amplitude 
sum distributions, poiQ ~), is simply obtained by 
transferring the output of the gated integrator 
directly to the MCA through amplifier A3 and gate 
G2. For this measurement, G2 is kept open contin- 
uously by proper positioning of S4. The amplifier 
A3 can be a combination of the amplifier built into 
the integrator and an external pulse amplifier. The 
gain of A3 is adjusted to give the desired range 
acceptable to the MCA, e.g., so that the maximum 
pulse amplitude is below 8 V. 

The distinction between the distributions 
/7i(^f |A<^*) and pi(qf\(}>f) in Table 4 is that 
A<^* is an arbitrary fixed phase window selected by 
DDGl whereas <f>f is a fixed phase associated with 
the occurrence of /th pulse. The measurement of 
Piiqf |<^f ) is like that iox piiqf |<^f ,j2') except 
that j2 * is not specified. In cases where specifi- 
cation of j2* is not required, the pulse selector 
circuit can be controlled directly by the out- 
put pulse from DDGl as for the case considered in 
Fig. 5. 

The measurement of/7i(^f jQ') is like that for 
piiqf |<^f ,<2*) except <^f is not specified. The 
measurement of j33(<^f |<A-i,9i-~i,<2*) is like that 
for p2((j>f I </»,•- i,^,-i) considered in Fig. 8 with an 
additional specification on the value for Q ', which 
is achieved by controlling the pulse selector circuit 
with DDG2 rather than DDGl. The measurement 
oipiiqf |A</>,*-i) is like that {oTpi(qi\Atj-i) in the 
constant excitation case with the exception that the 
At control logic is triggered by a particular pulse 
selected from the pulse selector circuit so that / has 
a specified value, i.e., it is not arbitrary as in the 
constant excitation case. For the measurements 
of the conditional distributions piiqf \<l>f ,qi-i) 
and pi^qf \(l)f ,(j)i-i,qiZ.i), the values of <^,-i and 
qi-i are restricted by the windows defined 
respectively by DDG4 and SCAl as for the 
measurement of J32(<^f |<A-i,?i-i) shown in 
Fig. 8. 

It is obvious that the system in Fig. 4 can be con- 
figured to measure other types of conditional 
distributions such as ;74(^f |<^f ,<ft-i,^/-~i.<2*)- 
However, the operation of the system has only 
been tested for distributions like those listed in 
Tables 4 and 5. The cycle selector allows triggering 
of DDGl and DDG4 only after n cycles have 
occurred. By using this feature it is possible to 
check for memory propagation from half-cycles 
that occurred prior to the most recent half-cycle 
such as would be indicated by the conditional 
distributions j3i(<^,^ \Ql) with k >j. 



4. Measurement System Components 

This section provides detailed information about 
the individual circuits that were designed for use 
with the measurement configurations shown in 
Figs. 3 and 4. Some of the circuits, such as the 
pulse-sorter circuit in Fig. 3 used for measurement 
ofpo(Atj) have previously been described [24] and 
remain unchanged. These circuits are not covered 
in this section. Some such as the At control logic 
circuits (parts A and B) have been revised and are 
included here. Qrcuits that have features specific 
to the measurement of phase-correlated distribu- 
tions such as the gated integrator, pulse selector, 
and gated single-channel analyzer in Fig. 4 have 
not previously been described and are covered 
here. Some circuits such as the zero-crossing detec- 
tor, pulse counter-cycle selector, and absolute- 
value selector are considered to have well known 
design and operating characteristics [31] and are 
not included here. Other circuits not considered 
here are those that are commercially available such 
as the digital-delay generators, time-to-amplitude 
converters, linear pulse amplifiers, and multi- 
channel analyzer. The description of each circuit 
given below includes both circuit and associated 
pulse diagrams. 

4.1 Time-Interval Control Logic (Parts A and B) 

The At control logic circuits described in this 
section have replaced the circuits previously shown 
in Figs. 2 and 4 of Ref. [24]. The designs have been 
improved to extend the measurement capabilities 
of the system and to eliminate or reduce errors 
previously noted [1]. 

4.1.1 Part A The function of the At control 
logic (part A) is to control a digital-delay generator 
(DDG3 in Figs. 3 and 4) to enable either the gate 
to the MCA (G2 in Figs. 3 and 4) or the gate for 
the stop pulse in the pulse-selector circuit of Fig. 9. 
The operation of the circuit can be understood 
with the aid of the pulse diagram shown in Fig. 10. 

Unlike the circuit previously described in Ref. 
[24], the buffer amplifiers defined by the transis- 
tors Tl and T2 are connected to separate inputs 
(f and e"). For some applications, such as the 
measurement of ;7i(^;|A/y-i), the inputs f and e" 
are connected together, and in other applications, 
such as the measurement of P2iqj\qj-i,Atj-i), 
these inputs are connected to different locations, 
e.g., by switch S2 and Fig. 3. The gain of these 
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Fig. 9. Diagram of the time interval (At) control-logic circuit, part A. The individual integrated circuit components are specified in the legend. 
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Fig. 10. Pulse timing diagram for the different indicated circuit 
locations in the A; control logic circuit (part A) shown in Fig. 9. 



buffers is adjusted so that their output is a 5 V 
logic pulse independent of the peak voltage of the 
input pulse. Ideally the input pulse voltage should 
lie within the range of 0.3 to 10 V corresponding to 
the range accepted by the MCA. For the systems 
shown in Figs. 3 and 4, the range of input pulse 
voltages is usually determined by the gain setting of 
the input amplifier Al. 

The outputs of the buffer amplifiers (Tl and T2 
in Fig. 9) trigger the 10-MHz serial shift registers 
SRI and SR2 that serve the purpose of delaying 
the pulses by the times n and t2 respectively. A 
pulse appearing at e" is allowed to trigger the 
digital-delay generator (DDG3) at output f" if the 
flip-flop circuit Fl defined by gates G5 and G6 is in 
the proper "initial-condition" state. If it is in this 
state, then the output of SR2 will cause it to 
change state after a delay of T2 = 600 ns, and the 
output of G5 will drop to zero thus triggering the 
one shot OSl. The one-shot produces a 2 /as pulse 
for triggering DDG3 and also sets flip-flop F2 
which indirectly results in the enabling of gate G12. 
If G12 is enabled, the returning pulse from DDG3 
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which appears at the f input will be allowed 
to pass to the output h that is used to control 
the gate to the MCA or the gate of the pulse 
selector. 

A pulse appearing at the f input will disable G12 
and prevent transfer of the DDG3 output pulse to 
h. In Fig. 10, the DDG3 pulse is indicated to have a 
delay At„-i and a width 5(Ar„-i). The circuit is 
thus designed not to record a pulse in the MCA if 
that pulse follows the pulse which triggered DDG3 
within a time less than A/«-i + T2. This feature 
ensures proper measurement of distributions con- 
ditioned on a fixed time separation, e.g., 
pi(qj\Atj-i), where it is required that A/j_i fall 
within the range Atj-t to Atj-t + 8{Atj-i) defined 
by DDG3. 

If a pulse appears at f, it will cause F2 to be reset 
after a delay of n = 300 ns set by SRI. The delay n 
must be less than T2 in order to prevent immediate 
resetting of F2 for coincident inputs at f and e". 
The pulse, generated by OS3, that initiates the dis- 
abling of gate G12 is delayed by a short interval 
defined by one-shot OS2. This small added delay 
(At' in Fig. 10) is necessary to ensure that the 
MCA gate stays enabled long enough to allow 
recording of any pulse that appears within the 
desired time interval. 

The circuit is reset by the pulse returned to input 
f by DDG3. The one-shot OS4 is triggered by the 
trailing edge of the DDG3 pulse independent of 
whether or not the gate G12 was disabled before 
time At„-i + S{At„-i) by an intermediate or 
recorded input pulse. The output of OS4 is delayed 
by the shift register SR3. Two of the outputs from 
SR3 control flip-flop F3 (gates G7 and G8) which 
in turn clears the contents of SR2 to prevent simul- 
taneous setting and resetting of Fl. Another out- 
put of SR3 then sets Fl to the initial condition 
which allows DDG3 to be triggered at output f by 
the next pulse to appear at e". Generally, any pulse 
that appears at e" will also appear at f so that F2 
will also be reset in the event that it was not 
already reset. 

The circuit has an additional, optional "gate" 
input which allows the output pulse at h to be 
gated by an independent source. This gate is also 
designated by G3 in Fig. 4 and used for the 
measurement of p2(qf \4>f yQ"^)- The present At 
control logic circuit has an auto-reset capability 
similar to that described previously [24] so that a 
reset pulse will appear at Gl to initialize Fl within 
a time of 0.8 s if for some reason a pulse is not 
returned from DDG3. If it is necessary to record 



times longer than 0.8 s, then the 10 Hz auto-reset 
clock can be replaced with one of lower frequency. 

4.1.2 Part B Part B of the At control logic is 
required together with part A when measurements 
of the conditional pulse-amplitude distributions 
pi{qj\Atj-i) are made for i>l as discussed in 
Sec. 3.1. The version of the circuit shown in Fig. 11 
differs from that described in our earlier work 
(Fig. 4 of Ref. [24]) which only enabled measure- 
ment of pi(^;|A/j_2) corresponding to the case 
where i=2. The present circuit incorporates a 
pulse counting feature that allows determination of 
distributions for i > 2. 

The operation of the circuit shown in Fig. 11 can 
be understood from a consideration of the pulse 
diagram shown in Fig. 12. The pulse time separa- 
tion between the (j —i + 1) and (j —i)th events (see 
Fig. 1) is restricted by the delay and pulse width 
settings of DDG3. The output of DDG3 is con- 
trolled by part A as described in the previous sec- 
tion so as to ensure that A/;_, is properly defined as 
a time separation between adjacent pulses. The op- 
eration of the At control logic (part B) is thus ini- 
tiated by the simultaneous occurrence of an event 
(PD) pulse at input eb = e" and a pulse from the 
h = ha (or equivalently SC2) output of the At con- 
trol logic (part A) at input SC2. An event occurring 
in coincidence with the SC2 input sets flip-flop F5 
(gates G14 and G15) which in turn allows the 
counter C02 to count subsequent event pulses. 
When the output of the counter corresponds to a 
binary number equal to i — 2, as determined by the 
settings of switches S1-S4, the comparator CI trig- 
gers the one-shot OSS. The output of OSS sets flip- 
flop F6 (gates G17 and G18) which enables the 
gate (G2 in Fig. 3) of the MCA at output h = hb. It 
also triggers the one-shot OS6 which in turn resets 
FS and consequently also the counter C02. The 
MCA will record the next (j th) event to occur after 
G3 is opened independent of its time separation 
from the preceding (j — l)th event. Immediately af- 
ter this event is recorded, the MCA sends a S V 
"event pulse" to the input d'. After a short delay 
determined by the one-shots OS7 and OS8, the 
event pulse resets F6 and thereby disables the 
MCA gate. The circuit is then ready to be triggered 
by the next coincidence to occur at inputs SC2 and 
eb = e". If for any reason the MCA fails to return an 
event pulse, e.g., the ; th pulse was too small to be 
recorded, then the counter C03 automatically pro- 
vides a reset for F6 after 8 successive coincident 
pulses have appeared at SC2 and eb = e". 
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Fig. 11. Diagram of the At control logic circuit, part B. The individual integratcd-circuit components are specified in the legend. 



1 


X 




n 


1 


2 


3 


4 5 


X 


X Event Pulse 




1 

1 
1 


"TV 












^f] 




UbDpul 












1 


^ 


<-AtbyOS7 


h 


\ 


II 


V MCA gate 


•^b 














1 
''' MCA event 










Count selector set at £ 







Fig. 12. Pulse timing diagram for the At control logic (part B) 
shown in Fig. 11. It is assumed in this example that the binaiy 
count selector is set by switches S1-S4 to the value 5. 



4.2 Pulse Selector 

A diagram of the pulse-selector circuit is shown 
in Fig. 13. The associated pulse diagram is shown 
in Fig. 14. This circuit is used to select a particular 



numbered event that occurs after a specified time 
or phase. It is used in the system shown in Fig. 4 to 
measure the conditional phase or amplitude distri- 
butions of specific pulses that occur within well- 
defined phase windows of the excitation voltage, 
e.g.,pi(<Af|e"),/ = 1,2,3... 

The pulse selector produces a - 2 V "start" pulse 
at output e that is coincident with the rising edge 
off a + 5 V digital-delay generator pulse applied to 
input c". This DDG pulse is assumed to define a 
fixed phase or time window, and therefore the 
output at e occurs at a known phase or time, i.e., it 
provides the appropriate phase reference point. 
The circuit is designed to produce - 2 V "stop" 
pulses at outputs e" (Stop 1) and SCI (Stop 2) that 
are coincident with the i th pulse to occur after the 
phase-reference point, i.e., after the "start" pulse. 
The value of i is determined by the setting of the 
binary count selector (switches SI - S8). 
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Fig. 14. Pulse timing diagram for the pulse-selector circuit 
shown in Fig. 13. It is assumed here that the binary count selec- 
tor is set to select the third pulse to occur within the time inter- 
val defined by the duration of the pulse appearing at input c". 



In the operation of this circuit, the pulse coun- 
ters COl and C02 are enabled by the DDG 
window. The negative event pulses that appear at 
input f are then counted by COl and C02. The 
binary outputs of these counters are sensed by the 
comparators CI and C2. If and when the binary 
number presented by the counter outputs equals 
the 8-bit binary number selected by the terminal 
switches, i.e., for i-\ to 256 pulses, the output of 
C2 goes from to -h 5 V and sets the flip-flop 
defined by gates Gl and G2. This flip-flop triggers 
the one-shot 0S3 that produces the "stop" pulses 
that ultimately appear at e" and SCI. 

At the falling edge of the DDG pulse, the coun- 
ters COl and C02 are reset and the one-shots OSl 
and 0S2 are triggered. The output of 0S2 resets 
the flip-flop and delivers a "stop" pulse to SCI in 
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coincidence with the end of the DDG window. This 
pulse does not appear at the e" output. The 
"end-of-window" stop pulse provides another 
phase mark that may be useful in some applica- 
tions. 

The appearance of the stop pulse at e" can also 
be controlled by another external pulse applied to 
the "gate" input. This option is used for the system 
shown in Fig. 4 for the measurement of various 
conditional phase distributions as indicated in 
Table 5. 

4.3 Single-Channel Analyzer 

The single-channel analyzer circuit designed for 
the measurement systems shown in Figs. 3 and 4 is 
presented in Fig. 15. The operation of this circuit is 
indicated by the pulse diagrams shown in Figs. 16 
and 17. If the amplitude of the input pulse at b' lies 
within a selectable voltage window, then the SCA 
delivers both a - 2 V and + 5 V pulse to the indi- 
cated output points at b". 



The circuit is capable of either gated or ungated 
operation. For ungated operation, output pulses 
are generated at b" for any input pulse that has an 
amplitude within the selected window independent 
of its time of occurrence. For gated operation, 
output pulses are generated only if the input pulse 
occurs within a time interval defined by the width 
of a 5 V gate pulse applied to input d". Depending 
on the setting of the switch S2, the output pulses 
for the gated operation will either appear at a time 
approximately coincident with the input pulse 
(with a slight delay) or at a time coincident with 
the end of the gate pulse. The three possible 
modes of operation are identified in Fig. 16. 

In the operation of this circuit, the event pulse is 
sensed by amplifier Al, the output of which either 
follows or inverts the signal depending on the posi- 
tion of switch SI. The output of this amplifier then 
proceeds to the analog comparators CI and C2. 
The other inputs to the comparators are derived 
from that part of the circuit (amplifiers A2-A5) 
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Fig. IS. Circuit diagram for the gated single-channel analyzer. The individual intcgrated-circuit components and switch configurations 
of the different modes of operation are shown in the legend. 
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possible operating modes for the single-channel analyzer shown 
in Fig. 15. 
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Fig. 17. Pulse timing diagram for the indicated circuit locations 
of the single-channel analyzer shown in Fig. 15. 



that define the minimum voltage and width of the 
window. Voltages in the range of to + 5 V are 
selected by the two 2kCl resistors denoted by 
"min" and "width" in Fig. 15. These voltages are 
doubled by amplifiers A2 and A3. Amplifier A4 
sums and inverts the outputs of A2 and A3, and 
amplifier A5 in turn inverts the output of A4. 
Consequently, the positive input of CI is a voltage 
between and + 10 V corresponding to twice the 
"min" value and the positive input of C2 is a 
voltage between and + 15 V equal to twice the 
value of the "min" plus "width" voltages. 
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As indicated in Fig. 17, the outputs of CI and C2 
are normally high ( + 6V) and go negative when 
the negative input exceeds the positive input. The 
circuit is designed so that input event pulses with 
amplitudes below the "min" value are ignored 
while those exceeding the maximum Value ("min" 
plus "window") are inhibited. In the latter case, 
the inhibition results from the setting of flip-flop 
Fl (gates Gl and G2) by the output of C2 which in 
turn disables G3. If the input event pulses fall 
within the window, then the output of the one-shot 
OSl triggered by CI passes through the gate G3 
and becomes the source for the 2 ju.s output pulse 
at b". The output of OSl also triggers OS2 which 
resets Fl after a delay sufficient to prevent passage 
of any pulses through G3 that exceed the maximum 
value as illustrated in Fig. 17. 

For normal gated operation, the switch S2 is set 
so that bo = bi. This allows G4 to be enabled by a 
+ 5 V gate pulse at d" and thereby permits passage 
of the pulse from G3 to the output buffer amplifiers 
(transistors T1-T3). If it is desired to have the out- 
put pulses appear at a fixed time corresponding to 
the end of the gate pulse, then bo=b2 at S2. In this 
mode, the output of OSl passes through G3 and G4 
and ultimately triggers flip-flop F2 (gates G5 and 
G6) which in turn enables 07. The one-shot 0S3 is 
triggered by the falling edge of the gate pulse and 
thus produces a pulse that passes through G7 and 
triggers OS4. The output of OS4 resets F2 and also 
becomes the source of the output pulses at b". 

4.4 Gated Integrator 

The integrator circuit used for the measurement 
system in Fig. 4 is shown in Fig. 18 and the corre- 
sponding pulse diagram is shown in Fig. 19. The 
output of the integrator is a pulse with an amplitude 
in the range of to 12 V directly proportional to the 
sum of the areas under all pulses that occur within 
the gate time interval denoted by A/i in Fig. 19. If all 
pulses have the same shape so that their amplitudes 
are proportional to their areas, then the height of 
the integrator output pulse is also proportional to 
the sum of the amplitudes of all pulses occurring 
within Ati. 

For the circuit in Fig. 18, the input pulses to the 
integrating amplifier Al are assumed to be of con- 
stant width ( ~ 1 ju.s) with amplitudes in the range of 
to - 12 V. In the absence of a 5 V gate pulse at 
input c, the flip-flop defined by gates Gl and G2 
keeps the field-effect transistors (FET's) Tl and T2 
turned on so that the 0.01 ju.F integrating capacitor 
is effectively shorted. The application of a pulse at 
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Fig. 18. Circuit diagram for the gated integrator. 
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Fig. 19. Pulse timing diagram for the gated integrator circuit 
shown in Fig. 18. 

c changes the state of the flip-flop which then turns 
off Tl and T2 thus allowing charge to accumulate 
on the integrating capacitor. At the end of the gate 
pulse, the one-shot OSl is triggered and its output 
momentarily turns on T3 which allows transfer of 
the integrator amplifier output voltage to amplifier 



A2 and thereby to the output terminal b'. The 
trailing edge of the OSl output pulse also triggers 
OS2 that generates a pulse to reset the flip-flop 
which in turn turns on Tl and T2 thereby discharg- 
ing the integrating capacitor. 



5. Examples of Results 

The purpose of this section is to show examples 
of data on conditional and unconditional distribu- 
tions that have been obtained using the systems 
shown in Figs. 3 and 4 respectively for measurement 
of dc and ac excited pulsating PD. A detailed 
discussion of the physical bases for the observed 
stochastic properties of PD phenomena goes 
beyond the scope of this paper and can be found in 
other works [1-4]. 
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5.1 Continuous Excitation Process 
(de-Generated PD) 

Shown in Fig. 20 are examples of the measured 
unconditional and conditional pulse-amplitude 
distributions poiqn), pi(?„ |Ar„-i), and 
/'2(?«|Af„-i,9„-i) for negative corona (Trichel) 
pulse discharges generated with a point-plane elec- 
trode gap in a neon-oxygen gas mixture at atmo- 
spheric pressure (100 kPa). The unconditional and 



first-order conditional distributions are plotted on 
a logarithmic scale and normalized to the maxi- 
mum values to facilitate comparisons of the various 
distributions. The pulse amplitudes are expressed 
in units of (pC) as explained in previous work 
[1,27] (also see Sect. 6.1). The dependence of the 
first-order distributions pi(qn\At„-i) on At„-i 
implies a strong positive dependence of q„ on 
At„-i, i.e., (At„-i T => ?n T)- This behavior can be 
explained in terms of the expected influence of the 
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Fig. 20. Measured unconditional and conditional pulse-amplitude distributions pa{q„), p\{q„\^t„-\), and 
P2((?n|i?n-i,A/„-i) at the indicated "fixed" values for A/„-i and q„.\ for negative-corona discharge pulses generated 
using a point-plane electrode gap in a Ne-^5% O2 gas mixture (see Ref. [1]). 
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moving negative-ion space-charge cloud from the 
previous pulse on the electric field in the gap and 
consequently also on the growth of the next 
discharge pulse [1]. 

The dependence of the distributions 
P2(qn\At„-i,q„-i) on^„-i for a fixed At„-i implies 
a negative dependence of q„ on (j„-i, i.e., 
(q„-i1 ^ q„ i,At„-i ). This behavior can be 
explained from consideration of the size of the 
space-charge cloud from the previous event on the 
growth of a discharge pulse. The data for both the 
first and second order pulse-amplitude distributions 
clearly demonstrate the importance of memory 
effects in determining the stochastic behavior of 
this discharge phenomenon. 

The dashed lines shown for the first-order distri- 
butions at A?„-i = 177, 197, and 217 |xs were calcu- 
lated using the integral expression [1] 



pi(q„\At„-i)=po(At„-i) ^ po(qn-i) 

Jo 



Pi(At„-i\q„-i)p2iqn\qn-i,Atn-i)(iq„-i (13) 

with numerical data obtained for the distributions 
shown in the right-hand side. Data for the condi- 
tional time-interval distribution,pi(A?„ \q„ ), used in 
the integral are shown in Fig. 21. It is interesting to 
note that in this case {At„(q„)) increases as q„ 
increases, i.e., (q„ t ^ At„ t). This means that the 
larger the previous event, the longer on average will 
be the time spacing between this event and the next 
event. This has been explained in terms of the 
influence of the electric field generated by space 
charge from earlier discharge pulses in suppressing 
the release of electrons from the cathode needed to 
initiate subsequent pulses [1]. 

It is evident from the results shown here that it 
would be impossible to find a physical interpreta- 
tion of measured unconditional pulse-amplitude 
distributions without information about the mem- 
ory effects revealed by the conditional distributions. 
The unconditional amplitude distribution is related 
to the time-interval distribution po(At„) and the 
conditional distribution pi(q„\At„-\) through 
Eq. (7). The first-order conditional distribution is in 
turn related to higher order distributions through 
Eq. (13). An unraveling of memory effects is a 
required step toward understanding the observed 
stochastic properties of random point processes 
such as reported here for the Trichel-pulse dis- 
charges. 



5.2 Periodic Time Varying Excitation Process 
(ac-generated PD) 

Data were obtained in this case for partial dis- 
charges generated by applying a sinusoidal alter- 
nating voltage to a point-dielectric discharge gap in 
air. Preliminary results from these measurements 
have recently been reported [2, 4, 32]. Figure 22 
shows examples of measured unconditional and 
conditional pulse-amplitude distributions of the 
first negative pulse to appear in each cycle. Also 
shown are the unconditional and conditional 
phase-of-occurrence distributions for this pulse. 
These results were acquired after observing the 
discharge pulses for many thousands of cycles of 
the applied voltage. 

The data shown in Fig. 22 were obtained using a 
stainless-steel point electrode positioned over a 
large, flat polytetrafluoroethylene (PTFE) dielec- 
tric surface in room air at a temperature of 23 °C. 
The tip of the stainless-steel electrode had a 
radius-of-curvature of 0.05 mm and was separated 
from the PTFE surface by a gap of 1.2 mm. A 
200 Hz, 3.0 kV rms sinusoidal voltage was applied 
to the gap. 

All distributions shown in Fig. 22 have been arbi- 
trarily normalized to the maximum values. The 
indicated values for Q * correspond to the integra- 
ted charge associated with all positive PD events in 
the previous half-cycle [see Eq. (4)] and define the 
type of line used to represent the data. In the case 
of the second-order distributions, /J2(?r| <^ 7,(2^ )j 
the fixed phase windows are defined directly under 
the data to which they apply. 

There are clear indications from these data of 
the significance of memory propagation in deter- 
mining the stochastic behavior of the phenomenon. 
The data for /Ji(<^r|j2'^) indicate that the larger 
the value of Q*, the smaller is the value of the 
mean phase-of-occurrence of the first negative PD 
pulse. This means that <f>T has a negative depen- 
dence on Q ■", i.e., (0 "" t => «^ I i). The data for 
P2(qT\<pi ,Q*) show that ^7 is positively depen- 
dent upon j2 ^ for a fixed phase-of-occurrence <f) T, 
i.e., (!2^t=>?rT,0r). The data for ;ji(^ri<^r) 
show that the mean value of the first negative pulse 
amplitude increases with its phase-of-occurrence. 
This distribution is related to the unconditional 
distribution po(Q*) and the other conditional 
distributions shown in Fig. 22 by the expression 



Po(<f>i)pi(qi |<^ 



■)= j Me") 

Jo 



Pii<}>T\QnP2iql\<t>T,Q^)6Q^ 



(14) 
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Fig. 21. Measured conditional time-separation distributionspi(A(„!9„) at the indicated values 
for Ai„ or q„ for negative-corona discharge pulses generated using a point-plane electrode gap 
in a Ne + 5% O2 gas mixture under conditions similar to those that yielded the data shown in 
Fig. 20 (see Ref. [1]). 



The corresponding data forpo(Q*) are not shown. 
At present, it has not been possible to obtain 
enough data on the required distributions under 
stationary discharge conditions to verify that 
Eq. (14) is indeed consistent with the experimental 
results. 



It has recently been shown [33] that the types of 
stochastic behavior for ac-generated PD reported 
here are consistent with theoretical predictions 
derived from a Monte-Carlo simulation of the phe- 
nomenon. The primary long-term (cycle-to-cycle) 
mechanism for memory propagation is that due to 
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electric charge accumulation on the dielectric 
surface during a PD event. It is well known that a 
quasi-permanent surface-charge distribution can 
exist on a solid insulating surface for times that are 
long compared to typical periods of the excitation 
voltage [34H37]. A significant fraction of the 
charge deposited on a dielectric surface by a PD 
event will remain to affect the local electric-field 
strength at the site where the next PD event is 
initiated. Both the probability for PD initiation and 
the distribution of the PD amplitudes depends at 
any given time on the local instantaneous electric- 
field strength. As in the case of dc-generated 
Trichel pulses, short-term pulse-to-pulse memory 
propagation can also exist for ac-generated PD. 
Mechanisms for memory propagation in this case 
could include moving ion space charge [1, 38], 
diffusion of metastable excited species [39], or a 
rapid redistribution of charge on a dielectric 
surface following a PD event [2, 40]. 

6. Calibrations and Sources of Error 

In the discussion about the earlier version of the 
stochastic analyzer [24], several sources of system- 
atic error were considered. These were primarily 
errors associated with the finite digital-delay gener- 
ator time window and the finite reset time of the 
time-to-amplitude converter. These among other 
sources of error need to be considered in making 
interpretations of the measured distributions and 
in judging the validities of consistency analyses 
performed using relationships like Eqs. (7), (8), 
and (13). It is, for example, important in consider- 
ing the use of Eq. (7) in checking consistency 
among the measured distributions po(?; ), po(Atj ), 
and pi(qj\Atj-i) to know the extent to which 
pi(^;|A^;-i) represents the true conditional distri- 
bution for a fixed Atj-i [see Eq. (9)]. It is also 
important that measurements of Atj using a TAC 
and the determination of A^;-i using the combined 
DDG and A^ control logic yield identical time 
separations. Any error in one of these circuits rela- 
tive to the other can cause difficulties in perform- 
ing the integration implied by Eq. (7). Thus, for 
example, it is generally necessary to make correc- 
tions for the delay T2 introduced by the A^ control 
logic circuit (see Fig. 10). In this section we 
consider the possible sources of error in the 
measurement of various amplitude, phase-of- 
occurrence, time-separation, and integrated pulse 
(charge) distributions that can be measured with 
the system described above. Methods for calibra- 
tion and testing of system performance are also 
discussed. 



6.1 Amplitude Distributions 

6.6.1 Pulse Shape Considerations The method 
for calibration of pulse amplitudes for PD has been 
described previously [27]. One could, in the 
simplest case, directly apply pulses of a known 
amplitude to the input of the system (amplifier Al 
of Figs. 3 and 4) and then record the MCA channel 
numbers corresponding to pulses of different 
amplitude. In most cases, however, it is desirable 
that the simulated input pulses used for calibration 
be similar in shape to those observed for the 
phenomenon of interest. This is especially required 
in the case of partial-discharge measurements 
where the amplitude of the recorded PD event is 
supposed to be proportional to the discharge inten- 
sity. It has been shown [1, 27] that, for the types of 
PD phenomena considered in the previous section, 
the recorded pulse amplitude is proportional to the 
net charge generated during the pulse provided the 
width of the impulse response for the detection 
system is very large compared to the intrinsic width 
of a typical discharge pulse. Under this condition, 
the shape of the recorded pulse is governed 
primarily by the impulse response of the detection 
circuit. The width of the impulse response for the 
detection system used to obtain the data in Figs. 
20-22 is approximately 1.5 n,s compared to a typical 
intrinsic PD pulse width of 1 to 11 ns. 

Pulses for some types of PD phenomena such as 
pulsating corona in air are known to have tails that 
are long compared to the 1.5 \x.s impulse response 
width corresponding to the conditions under which 
the data in the previous section were taken. In such 
cases, not only is the measured pulse amplitude no 
longer directly proportional to the total charge 
generated by the PD event, but there may also exist 
the possibility that the system will sample the tail 
of the pulse one or more times in addition to its 
peak value. As noted in our earlier work [27], this 
problem can occur if the system sampling rate, 
governed primarily by the dead time of the MCA, 
is sufficiently high. (The MCA dead time is approx- 
imately 2 (AS for the system used in this work). If 
this problem occurs, the measured amplitude 
distribution will be artificially enhanced at the low- 
amplitude end due to recording of the tails. A 
similar problem of pulse amplitude distribution 
distortion is known to occur in cases where the PD 
occurs as bursts of pulses in which: 1) the spacing 
between pulses is comparable to or shorter than 
the detector impulse response time or MCA dead- 
time, i.e., the system sampling time; 2) the duration 
of the burst is comparable to or longer than the 
system sampling time; and 3) there is a high degree 
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of correlation among the amplitudes of pulses 
within a burst. Such short-duration, burst-type PD 
pulses for which pulse amplitudes are highly corre- 
lated are known to occur under some conditions 
[4, 27, 41]. 

The accidental sampling of pulse tails can also 
be a problem when measuring conditional pulse- 
amplitude distributions as discussed in our earlier 
work [24]. In this case, a problem arises if the 
MCA is gated on by the At control logic at pre- 
cisely the time when an event pulse is decaying, i.e., 
after a peak has occurred. The problem has been 
minimized in the present system by a combination 
of shaping the pulses that enter the MCA and by 
minimizing the delay time A^ indicated in Fig. 10 
as determined by the one shots 0S2 and 0S3 in 
Fig. 9. These adjustments were sufficient to yield 
acceptable results for the types of PD phenomena 
considered in this work. Complete elimination of 
this problem is difficult, but can be at least partially 
achieved by modification of the At control logic so 
that it senses if a pulse has occurred within a short 
time of approximately one pulse width before the 
DDG3 pulse is returned to input f" in Fig. 9. If a 
pulse does occur within that time, then a condition 
can be set to force a delay in the opening of the 
MCA gate at output h. 

6.1.2 Effect of Amplifier Nonlinearities Another 
possible source of distortions in measured condi- 
tional or unconditional pulse-amplitude distribu- 
tions is that associated with nonlinearities in the 
gains of pulse amplifiers used in the detection cir- 
cuitry or elsewhere in the measurement system, e.g., 
amplifier Al in Figs. 3 and 4. Figure 23 shows typi- 
cal examples of calibration curves used in analyzing 
the data on pulse-amplitude distributions such as 
those shown in Figs. 20 and 22. Shown are plots of 
amplitude in charge (pC) versus MCA channel 
number for two different ranges of amplitude and 
for two different gain settings (^i and ^2) of the 
amplifier Al for the same input amplitude range. 
The onset of gain saturation in the preamplifier that 
detects a pulses is indicated by the vertical arrow 
pointing to a place on the curve corresponding to 
the highest gain and highest amplitude range. 

It is desirable that the response of the pulse 
amplifier be as linear as possible. Deviations from 
linearity introduce complications in determining 
the true amplitude distribution from the measured 
data as will be shown below. The data recorded by 
the MCA can be represented by the array of 
numbers N(k) where k is the channel number and 
N(k) is the number of events stored in k. For a 256 



200- 



150- 




Fig. 23. Examples of amplitude-calibration curves for two 
different ranges of amplitudes and two different amplifier gains, 
gi, and g2, that apply to the measurement of PD pulse-ampli- 
tude distributions. The vertical arrow indicates the onset of gain 
saturation for the highest range and gain^i. 



channel MCA, k is restricted to integer values 
between 1 and 256. From calibration data such as 
shown in Fig. 23, one can relate the N(k) data to 
the true amplitude distribution by 



Pj(q)dq (15) 

a(ki ) 



where A^o is a normalization constant and pj(q) 
denotes a "true" / th order conditional amplitude 
distribution. If a single channel width is assumed so 
that Ak = l, then q (ku) and q (A;l) are defined here 
to correspond respectively to the values ku=k+^ 
and kL=k-\. 
If the pulse amplification is linear, we can write 



k = ao + aiq. 



(16) 
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where ao and ai are constants. If /?;(?) is slowly 
varying over the interval [q {ku), q (A:l)] so that 



and 



;>;(g(M)-P;(g(^L)) 
Pi (a) 



<l 



(17) 



where q =[q(Jcv)+q{ki)'\l2, then in the linear 
case 



Nik) « NoPi (q )[q (kv) -q^k^] (18) 



= ^o(^>,(^), 



(19) 



Equations (18) and (19) imply that N{k) versus k is 
a discretized approximation to the true distribution 

Pj(q)- 

For cases where the response is nonlinear, 
Eq. (19) is not valid and even Eq. (18) may fail to 
be a good approximation. If a quadratic depen- 
dence is included in Eq. (16) by adding a term aiq^ 
to the right-hand side, then the factor of (1/ai) 
in Eq. (19) must be replaced with the factor 
[l/(ai + 2a2'q )] which depends on ^. For a 
sufficiently large quadratic contribution, it is 
necessary to consider this "q dependent factor in 
attempts to estimate the;7;(^) distribution from 
the raw N(k) versus k data. 

Under conditions of severe nonlinearity, it may 
become impossible to extract meaningful informa- 
tion about pj (q ) from the MCA data. One such 
case is that encountered when the input amplifier 
gain approaches saturation. As an example of this 
case we consider an amplifier that begins to satu- 
rate for q>qs- The calibration curve for this case is 
represented mathematically by: 

k=ao + aiq,q<qs (20) 



k = ao + ai[l-exp{-pq)],q>qs, (21) 



where aj and /3 are constants and the gain is 
assumed to be linear for ^ <^s. In order that both k 
and dkldq be continuous at q =qs, the coefficients 
must satisfy the relationships 



ai = 



aiqs 



[l-exp(-pq,)] 



(22) 



(Pqs+l)exp(-Pqs) = l. 



In the region of saturation (q ^q^) we have 

dk 

-^ = aiexp[-p(q-q,)]. 



(23) 



(24) 



Equation (24) implies that for a fixed increment of 
channel number, e.g., AA: = 1, the difference 
between q (ku) and q (ki) increases exponentially 
with q, i.e.. 



q(ku)-q{kL)~exp[P{q -qs)]. 



(25) 



In this case, Eq. (18) may not hold and, in general, 
one must resort to Eq. (15) to relate N(k) to pj(q). 
The range of q values over which Eq. (15) must be 
integrated can become large enough under some 
conditions to prevent the determination of reason- 
able estimates for pj(q) from the MCA data. 
Generally, the effect of amplifier saturation is to 
cause N(k) to become artificially enhanced at large 
values of A:. The problems associated with satura- 
tion can usually be avoided by making careful 
adjustments of amplifier gain. 

A reasonable estimate of the true amplitude 
distribution also requires that the amplitude incre- 
ment, q (ku) —q (ki), associated with the width of a 
single channel be small compared to the character- 
istic width of the distribution. This can usually be 
assured by appropriate adjustment of the amplifier 
gain and MCA pulse discrimination level that re- 
spectively determine the parameters ai and ao in 
Eq. (16). There may exist cases, however, where 
the ability to make a precise determination of con- 
ditional or unconditional distributions is severely 
limited by the inherent resolution of the MCA. 

6.1.3 Noise Broadening Under low signal-level 
conditions, distortion of the measured amplitude 
distributions can result from effects of noise. 
Specific sources of noise are not identified here, 
but they could simply be those associated with 
normal amplifier operation. The type of noise 
considered in the present discussion is often 
referred to as "white noise." The noise is assumed 
to have constant statistical characteristics during 
the time of a typical measurement, i.e., it is as- 
sumed to be stationary. 

Not considered in this discussion are erratic or 
time-dependent noise such as might appear as ran- 
dom or phase correlated pulses generated by 
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pick-up from sources external to the system of 
interest. In cases where such externally produced 
impulses cannot be eliminated by adjusting the 
amplifier or MCA discrimination levels, it is possi- 
ble that these impulses will introduce severe distor- 
tions, especially if they are narrowly distributed in 
amplitude, phase, or frequency. Elimination of 
interference from impulse noise sources can be 
achieved under some conditions by using shielding 
or digital-filtering techniques [42, 43]. Discussion 
of these techniques goes beyond the scope of the 
present work. 

The imposition of a constant background noise 
on the detected impulse signals introduces a broad- 
ening in the amplitude distributions recorded by 
the MCA. The broadening effect can be estimated 
from the convolution 



Piiq 



>\: 



Pjiqlfiq -q')<iq'. 



(26) 



where Pj{^) is the broadened distribution and 
fCq —q') is a function that represents the statisti- 
cal distribution about a mean value, 'q , due to 
noise. In cases where the noise is inherent in the 
measurement system, the form of f{q —q') can 
sometimes be estimated from the calibration data. 
Figure 24 shows an example of a set of calibration 
data recorded by the MCA under conditions where 
pulses of known amplitude, ^,(/ = l,2,...), are 
injected from the calibration source during fixed 
intervals of time. If noise were not present, counts 
would be recorded in only one channel of the MCA 
for each value of qi . The fact that counts appear in 
12 or more channels for the data shown in 
Fig. 24 means that there is some broadening due to 
the presence of noise. In most cases like that 
shown in Fig. 24, the noise can be approximated by 
a Gaussian function, i.e., 

fiq -q') = i'rrwy'cxp[-{q -qy/w], (27) 

where the width, w, is independent of "q . Broaden- 
ing due to noise can be significant if the condition 
A^ >w is not satisfied, where A^ is the characteris- 
tic width of pj(q). If noise broadening is deter- 
mined to be significant, then it may be possible to 
develop a deconvolution procedure using Eqs. (26) 
and (27) to obtain an improved estimate of pj{q) 
from the measured data recorded by the MCA. No 
attempts have been made to implement noise 
deconvolution procedures in the present work. 
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Fig. 24. Amplitude-calibration data recorded by the MCA that 
show broadening due to noise. 



6.2 Amplitude Sum (Integrated-Charge) 
Distribution 

The measurement of integrated-charge distribu- 
tion, /io(!3 ")) for a specified phase region (positive 
or negative half-cycle) is subject to the same errors 
considered above for pulse-amplitude distribution, 
e.g., effects of amplifier nonlinearities and noise 
broadening. In addition to these, there are other 
mechanisms for introducing systematic errors that 
can result from the operating characteristics of the 
gated integrator and its associated input amplifier. 
These possible sources of error are examined here. 
Methods for calibration of the integrator that can 
reveal systematic errors and precautions that can 
be taken to ensure proper operation of this circuit 
are also considered. 

It should first be realized that the output of the 
pulse amplifier (A2 in Fig. 4) is a rectangular pulse 
of constant width independent of the shape of the 
input pulse. The amplitude of the output pulse is 
directly proportional to the amplitude of the input 
pulse provided the input pulse lies above a critical 
value qj (discrimination level). The output of the 
amplifier, q, is related to the input, q ', by 



q =gA2q',q'>qd 
q = 0,q'<qi, 



(28) 
(29) 



where gA2 is a constant corresponding to the gain of 
A2. 

The input amplifier A2 essentially acts like a 
peak detector. By using this type of amplifier as an 
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input to the integrator, the output pulse of the 
integrator is forced to be proportional in amplitude 
to the sum of the amplitudes of the pulses appear- 
ing at the input to A2 consistent with Eq. (4). The 
integrator, therefore, does not yield an output that 
is a true measure of the integrated current associ- 
ated with the event pulses as given by 



G*=2 n'ifindt'. 



(30) 



where If(t') is the instantaneous current of the 
j th pulse and tj is the duration of the pulse defined 
such that If (t')^O only for times in the interval 
t + Tj>t'>t-Tj. The value of j2* is directly 
proportional to that given by Eq. (4) only if 
all pulses have the same shape. The error in the 
determination of true integrated charge is approxi- 
mately given by the difference 



^Q^=l\qf- n'lnndt'l (31) 



where the values for ^ fare determined by the cali- 
bration of amplitude in terms of charge-per-pulse 
under conditions where pulse shape is governed 
by the detector impulse response (see previous 
section). 

For cases where there may be events having 
pulse durations, t,-, that exceed the width of the 
impulse response, the sign of the error AQ- is 
most likely negative, i.e., the measurement prefer- 
entially tends to underestimate the integrated 
charge. An additional contribution to a negative 
error can occur if there are events with amplitudes 
lower than the discrimination level qa defined in 
Eq. (28). The integrated charge for these events 
is simply not included in the sum of amplitudes, 
Eq. (4). 

Calibration of the integrator requires the use of 
a gated pulse generator that produces a burst of 
pulses only during the time that the integrator is 
gated on as shown in Fig. 25. If Ni calibration 
pulses of known amplitude qc are applied during 
the gate interval, then the output pulse of the inte- 
grator should ideally have an amplitude directly 
proportional to the product Niqc. This means that 
the output should be independent of the number of 
pulses used if qc is maintained at a value of VdNj 



where Vc is constant. Thus, if Ni is increased by a 
factor of two, the integrator should give the same 
output provided q^ is correspondingly reduced by a 
factor of one half. 




Calibration Pulses 



Q = 5q 



Integrator Output 

Fig. 25. Pulse timing diagram that applies to a calibration of ihc 
gated integrator. The output of the integrator should ideally be 
a pulse with an amplitude directly proportional to the sum of 
the amplitudes of the calibration pulses that occur during the 
gale pulse. 

To ensure proper operation of the integrator, 
the dc bias level at the output of A2 should be 
adjusted to zero. After this adjustment is made, the 
offset voltage of the integrator (amplifier Al in 
Fig. 18) should be set at a value that forces the 
integrator output pulse amplitude to be zero when 
no pulses are applied to the input, i.e, when M = 0. 
This latter adjustment is required because the 
0.01 (Ji-F integrating capacitor in Fig. 18 can acquire 
an initial small charge attributable to the transient 
voltage associated with the opening of the FET's 
Tl and T2. Failure to make these adjustments will 
allow systematic errors to occur in the measured 
integrated charge due to an offset, i.e., a nonzero 
intercept of the calibration curve. 

Errors introduced by the presence of a finite 
offset are illustrated by the examples of calibration 
data shown in Figs. 26 and 27. Plotted in Fig. 26a is 
the amplitude (or equivalently charge) per calibra- 
tion pulse, qc = Q~/N], versus the recorded MCA 
channel number per pulse, N(k), for A^i = l to 5. 
Figure 26 shows the corresponding calibration plot 
of Q~ versus N(k), again for N] = l to 5. The data 
in Fig. 26a tend to fall on a straight line given by 



QVNi = 110+ ViN(k)/Nu 



(32) 



where tji is the slope and ijo is the intercept 
corresponding to a constant offset voltage. The 
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corresponding integrator calibration curves are 
given by 



Q-=Ni7]o + ViN{k). 



(33) 



500 




100 150 

N(k) / N, 



200 



250 



Fig. 26. Example of an integrator calibration performed under 
conditions where there is a finite offset corresponding to a non- 
zero intercept tju as shown in a). The calibration data in b) 
indicate the resulting uncertainty due to this offset (range of 
Q * defined by the solid lines). 



Consistent with the data shown in Fig. 26b, the 
intercepts for the (2* versus N(k) curves increase 
with increasing A^i. Under these operating condi- 
tions, there is an uncertainty of T}oi^Ni in the 
amph'tude sum (integrated charge) due to the off- 
set, where ANi is given by 



iiNi=Ni{max)-Ni{min) . 



(34) 



Here, Ari(max) and Ari(min) are respectively the 
maximum and minimum number of events that are 
likely to be recorded within the integrator gate 
interval (Afi in Fig. 19). The solid lines in Fig. 26a 
represent the error limits for the case considered 
when Ari(max) = 5 and Ni{mm) = 1. Figure 27 shows 
how this systematic error can be reduced if care is 
taken to minimize the offset so that 770=0.0. In 
general, all stages of amplification should be 
adjusted individually to eliminate offsets. Further 
reduction in the error can be achieved by giving the 
greatest weight to calibrations made using values for 
Ni that equal the mean number of experimentally 
observed events within the integrating interval. 
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Fig. 27. Example of an integrator calibration performed under 
conditions of zero offset (zero intercept, tjd) as shown in a); and 
corresponding error reduction as seen in b). 



6.3 Time-Interval and Phase-of-Occurrence 
Distributions 

Systematic errors that can occur in the measure- 
ment of pulse time-separation distributions were 
analyzed in our earlier work [24]. It was noted that 
the measured time-separation distributions can 
become distorted if: a) a significant fraction of the 
time separations are less than the time-to-ampli- 
tude converter reset time, Afr, and b) there are cor- 
relations among successive time separations. The 
reason why distortions are introduced under these 
conditions can be understood from a consideration 
of the example illustrated in Fig. 28. It is assumed 
that the phenomenon of interest appears in the 
form of pulse bursts where there is an ordering of 
pulse time separations within a burst such that the 
first separation is on the average smaller than the 
second and so on. This type of behavior occurs, for 
example, in the case of burst-type positive-corona 
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pulses generated in sulfur hexafluoride using point- 
plane electrode gaps [27]. 




AtiAt2At3 At4 Atg 



At,, 




-Atr-*^ 



b) 




Atr 



Fig. 28. Possible distortion of a measured time-separation dis- 
tribution, po(At„ ), for pulse bursts due to the finite TAG reset 
time Atr. 



Figure 28a shows pulse diagrams for two dif- 
ferent bursts of 7 pulses each with the indicated 
successive time separations Ari, Ar2,... and the 
TAC reset time Afr. It is assumed that two TACs 
are used for the measurement as shown in Fig. 3. 
The shaded time separations are those actually 
recorded by the system, i.e., Ati, and Ats are 
measured by TACl and At2 and Afj are measured 
by TAC2. The time separations Afa and AIa are not 
recorded because they occur before the TAC has 
had time to reset. The failure to record these 
separations causes the measured distribution to 
deviate from the true distribution as shown in 
Fig. 28b. This limitation can be overcome to some 
extent by using multiple TACs as discussed in the 
next section. 

In the case of conditional time-interval or phase- 
of-occurrence distributions, distortions can also 
occur if the range of values for the "fixed" variable 
are not sufficiently well restricted. As noted above 



[see Eq. (9)], the necessity of using a finite window 
size for the fixed variable introduces a broadening 
of the distribution. The problem has already been 
noted for the measurement of conditional pulse- 
amplitude distributions [24]. Unlike the broaden- 
ing due to noise, the broadening introduced by a 
finite window can be as3anmetrical with a resultant 
apparent shift in the associated mean value. Exam- 
ples of asymmetric broadening due to an increase 
in window for the variable Q ^ are shown in Fig. 29 
for measured conditional distributions /7i((^r|2"^) 
corresponding to the phase of the first negative PD 
pulse generated in a point-dielectric discharge gap. 
In general, it is desirable to keep the window size 
to the minin)um required to obtain acceptable 
statistics within a reasonable time. Excessive dis- 
tortions due to finite window size will invalidate 
consistency analysis using the various integral 
relationships among measured conditional and 
unconditional distributions, e.g., using Eqs. (7), (8), 
(13), and (14). 

7. Limitations, Extensions, and 

Alternatives 
7.1 Limitations 

It was noted at the outset that the system 
described in this work is optimally designed to 
investigate the stochastic properties of ac or dc 
generated PD pulses with repetition rates between 
50 and 5 x 10''/s. The lower limit on pulse rate is 
determined by the acceptable times within which 
observations can be made that will yield statistically 
significant results. If the phenomenon of interest is 
stationary over the time of observation, then there 
is, in principle, no lower limit on the pulse rate that 
could be observed. However, if the time between 
pulses exceeds the range of a TAC (typically 1 s), 
then time intervals must be measured by another 
method, e.g., using a digital clock with a gated 
pulse counter. For extremely low pulse rates (less 
than 1 per min) it may be more efficient to simply 
perform a statistical analysis of recorded data. 

The upper limit on the rate of observed pulses is 
imposed by time restrictions inherent to the elec- 
tronics. Specifically, there are limits due to: 1) the 
finite reset time of the TAC, 2) the dead-time of 
the MCA, 3) built-in delays in the SCA and At 
control logic circuits, and 4) the inherent impulse 
response of the pulse detection and amplifier 
network. It may be possible by using a broad-band 
detector and faster electronics to increase the 
range of applicability by an order-of-magnitude, 
i.e., to 5xl0'/s. However, this may introduce 
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Fig. 29. Examples of data for measured conditional phase-of-occurrence distributions 
Pi.(<I>T\Q*) that show asymmetric broadening due to a finite window for the "fixed" variable 



added complexity and cost that would make alter- 
native methods appear more attractive such as the 
"software" approach considered in Sec. 7.3. 

Although it was assumed here that pulse ampli- 
tude (or the sum of successive pulse amplitudes) 
is the appropriate "mark" for characterizing the 
intensity of the phenomenon, there may exist 



cases where other marks such as pulse area or 
pulse-shape parameters are more appropriate indi- 
cators of "intensity." The system documented here 
may still be applicable to these cases provided the 
mark can be converted to a pulse with an ampli- 
tude that is directly proportional to the "size" of 
the mark. 
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In principle, there is no upper limit to the size of 
a mark that can be measured. In some cases it may 
be necessary to restrict the amplitude of an event 
pulse by a linear attenuation network so that it 
does not exceed the voltage range acceptable to the 
MCA (0 to 8V in the present system). If the 
dynamic range of pulse amplitudes is very large 
(two or more orders-of-magnitude), then it may be 
necessary or desirable to replace the linear input 
pulse amplifier with a logarithmic amplifier. 

The lower limit on acceptable pulse amplitude is 
simply governed by signal-to-noise ratio. As the 
broadening due to noise (see Sec. 6.1.3) becomes 
comparable to the width of the observed distribu- 
tion, it becomes increasingly difficult to extract 
meaningful information about memory effects from 
the data. 

It was noted previously that when broadening 
due to noise or restricted variable window size is 
significant, it may no longer be possible to perform 
a consistency analysis among various measured 
distributions using the integral relationships that 
connect these distributions, e.g., Eq. (7). It may 
still be possible, nevertheless, to use the data on 
conditional distributions to establish the existence 
of memory propagation. The existence of memory 
propagation can be unequivocally established if it 
can be shown, for example, that the conditional 
distributions pi{qj\dktj-i€.{Ata, Atb)) and 
pi{qj\A.tj-i e (Ata; Atb')) do not coincide under at 
least one condition where the corresponding 
ranges of the time intervals (A/a, Atb) and 
(At„', Atw) are different. In order to determine this 
lack of coincidence, it is necessary: 1) to acquire 
enough data to demonstrate a statistically signifi- 
cant difference between the two distributions, and 
2) to acquire the data under conditions where the 
phenomenon is stationary. The effects of non- 
stationary behavior can be minimized if the data 
for the two conditional distributions can be 
accumulated simultaneously. Unfortunately, in the 
present system, which has only one At control logic 
circuit and one MCA, it is not possible to make 
simultaneous measurements of two distributions of 
the same type. The system can be operated, 
however, to alternately accumulate data in two 
different 256-channel segments of the 1024- 
channel MCA for two different ranges of the fixed 
variable. By periodic switching back and forth 
between the two segments, it may be possible to 
"average out" effects of nonstationary behavior. 
The existence of nonstationary behavior can often 
be detected from periodic monitoring of uncon- 
ditional distributions such as/7o(9/) or po{Q~) for 



which statistically significant data can be accumu- 
lated in much shorter times than for conditional 
distributions for which pulse count rates are lower. 
It has been shown [1, 4] that the profiles of uncon- 
ditional distributions tend to be more sensitive to 
nonstationary behavior than the profiles of condi- 
tional distributions. 

It should be realized that, in general, the time 
required to obtain statistically significant data for a 
distribution can increase rapidly as the number of 
restrictions imposed by the fixed variables is 
increased and as their ranges are reduced. This is 
perhaps the most stringent limitation inherent not 
only to this measurement system but also to 
stochastic analysis in general. For the types of 
pulsating PD phenomena investigated with this sys- 
tem, it has usually not been possible to obtain 
enough data with adequate statistics for con- 
ditional distributions higher than second order. 

7.2 Extensions 

Some of the limitations mentioned in the previ- 
ous and earlier sections can be overcome (or at 
least reduced) by introducing various extensions or 
expansions to the existing measurement system. 
One obvious extension would be to introduce addi- 
tional MCAs with associated circuitry to allow 
simultaneous measurement of two or more distri- 
butions. This would not only reduce the total data 
acquisition time, but would also allow better moni- 
toring of effects due to nonstationary behavior. 
Since nonstationary behavior in PD phenomena is 
often a consequence of discharge-induced modifi- 
cations of the discharge gap, e.g., changes in rates 
of electron emission from surfaces, multiple MCA 
measurement capability might make the system 
more useful as a diagnostic of insulation aging. 

The previously noted limitations on the measure- 
ment of time or phase separation distributions im- 
posed by the finite TAC reset time can be 
overcome by incorporating more TACs in series 
with an associated pulse sorting circuit so that each 
TAC measures a different successive time interval. 
If instead of using only two TACs one uses 2n 
TACs, where n>2, then the minimum time separa- 
tion that can be measured without introducing 
errors (see Sec. 6.3) is reduced from At,/2 to 
At,/2n. 

Of course, the introduction of added instrumen- 
tation such as MCAs and TACs significantly 
increases the cost of the system. In some cases this 
added cost may be more than compensated for 
by the reduction in time required to acquire and 
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analyze the data. Additional time savings may be 
achieved by automating the system to allow not 
only the simultaneous measurement of more than 
one distribution but also real-time calibration, data 
analysis, and optimization of time spent where 
needed to obtain the best statistics. 

Finally, it should be noted that although the 
system described here has only been applied to the 
measurement of pulse sequences in real time, it is 
also possible to use it for stochastic analysis of pre- 
recorded pulses. This merely requires that the 
source of the input pulses be derived from an elec- 
tronic recording device operated in the "play-back" 
mode. As will be argued in the next section, it 
should be possible to use prerecorded data derived 
from computer simulations to test the overall sys- 
tem performance. It is, of course, desirable that the 
simulation produce pulses with known stochastic 
properties, i.e., conditional distributions that mimic 
those of the phenomena under investigation, and 
that simulated pulses have an amplitude and shape 
similar to experimentally observed pulses. The use 
of the system described here to perform analysis on 
prerecorded data allows an obvious extension to 
pulses with repetition rates that are higher or lower 
than the ranges which are normally acceptable. 
However, one should consider whether or not it 
may be more efficient to analyze prerecorded 
data directly using computer software such as con- 
sidered in the next section. 

7.3 Alternatives 

One of the primary advantages of the system 
described in this work is the ability to measure the 
stochastic properties of a pulsating phenomenon in 
real time. It allows one to view on a computer out- 
put device (video monitor) the development of 
conditional or unconditional distributions as the 
data are acquired. With this capability, it is possi- 
ble for the operator to determine quickly the exis- 
tence of memory effects and to make decisions on 
the conditions that should be selected to yield the 
most interesting data. 

The system, in its present form, is designed to be 
a research tool for use in investigating memory 
propagation in pulsating phenomena. It is assumed 
that it is operated by those who have a thorough 
understanding of the phenomenon under investiga- 
tion. Because this is a highly "interactive" system, 
the quality of the information acquired from the 
measurements will be determined to a large extent 
from the judgments of the operator. It may be 
possible to construct an "automated" real-time 



stochastic analyzer based on the measurement con- 
cepts introduced here. Some aspects of the present 
system could, for example, be incorporated into 
advanced partial-discharge measurement systems 
that would allow the possibihty for meaningful 
pattern recognition needed to identify types of 
discharge sites. 

Another advantage of a real-time measurement 
system is that it overcomes problems of storing 
large data files. For example, the measurement of 
some of the second-order distributions shown in 
Figs. 22 and 23 required ten or more minutes of 
data acquisition time. This means that only a very 
small fraction of the total number of discharge 
events that occurred during the measurement time 
were actually recorded. A record of all discharge 
events that occurred in a typical 10 min segment 
would generally contain data for more than 10' 
pairs of numbers. Within the times required to 
obtain reasonably good statistics on higher-order 
conditional distributions, it is possible to generate 
data files containing all events that exceed mini- 
computer storage capacity. Analysis could then 
only be performed using either "main-frame" type 
computers or appropriately segmented data files in 
smaller computers. 

The obvious disadvantage of the present mea- 
surement system is that it does not make efficient 
use of the available data. In the measurement of 
conditional distributions, most of the information 
about the impulse events is discarded. Once the 
data are discarded, they can no longer be retrieved 
for subsequent analysis. The expense of introduc- 
ing additional MCAs in parallel to enhance the 
information retrieval efficiency can make the cost 
of the system prohibitive. 

In cases where either the amount of data is 
severely limited or the phenomenon is highly non- 
stationary, it may be essential to consider all of the 
available data. There may also be other cases 
where it is necessary to work with prerecorded data 
due to externally imposed geometrical or time 
constraints. 

The most efficient use of available data in such 
cases can, at least in principle, be achieved using 
an alternate approach that places more reliance on 
computer software. An example is given below of 
an algorithm that was developed to acquire data on 
conditional and unconditional distributions from a 
Monte Carlo simulation of ac-generated partial 
discharges. Details of the theoretical model upon 
which the simulation is based are given elsewhere 
[33] and will not be covered here. It need only be 
said that the simulation produces a sequence of 
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phase-correlated pulses with stochastic properties 
similar to those observed for PD generated by 
applying an alternating voltage to a point-dielectric 
electrode configuration [3]. 

A partial listing of a FORTRAN-77 routine used 
to sort data in real-time for determination of the 
distributions that apply to simulated PD is given in 
Table 6. The specific distributions considered in 
this table are the unconditional and conditional 
phase-of-occurrence distributions for the / th nega- 
tive PD pulse, i.e.,p[){(f>T) andpi{(f>T \Q*). 

Table 6. Partial listing of a FORTRAN 77 routine for deter- 
mining the distributions paC"^ 7 ) and/ii(i^ 7 |Q *) from numeri- 
cal data on 1^ 7 and two windows for Q *, i.e., Q* £ (Qitj, Qit) 
andQ*e(Q2t,Q2t). 



C POPI = ARRAY FOR UNCONDITIONAL PHASE DIST. 

C PIPII,PIPI2 = ARRAYS FOR CONDITIONAL PHASE DIST. 

C QPOS = INTEGRATED POSITIVE CHARGE (AMPLITUDE SUM). 

C QPUl , QPLI = LIMITS ON WINDOW-1 FOR QPOS. 

C QPU2 , QPL2 = LIMITS ON WINDOW-2 FOR QPOS. 

C IN = NEGATIVE PULSE INDEX NUMBER. 

C CHANNEL = INTEGER BETWEEN AND 200. 

C PHASE = PHASE OF NEGATIVE PULSE. 

CHANNEL = NINT(200'PHASE) 

POPI(CHANNEL,IN,0) = POPI(CHANNEL , IN , 0) + I 
IF(QPOS . LT . QPUI) THEN 

IF(QPOS.GT.QPLl) THEN 
P1PII(CHANNEL,IN,,0) = P1PII(CHANNEL, IN ,0) + I 

END IF 
END IF 
IF(QPOS . LT . QPU2) THEN 

IF(QPOS . GT . QPL2) THEN 
PIPI2( CHANNEL . IN , 0) = PIPI2(CHANNEL, IN.O) + 1 

END IF 
END IF 



The first statement in this routine converts the 
value of the normalized phase {4>1 /Inr) to the 
nearest integer value between and 200. It thus 
performs essentially the same discretization of the 
data as in a 200-channel MCA. The assigned 
integer value is then used to identify elements of 
three-dimensional integer arrays corresponding to 
particular distributions. In the case of the uncondi- 
tional phase distributions, the array element is 
increased by 1, i.e., by one pulse count. In the case 
of the conditional distributions, the array elements 
associated with pi{4>7 \Qu) and pi{(f>7 \Ql) are 
incremented by 1 only if the values of Qv and Q t 
lie within specified ranges. An example of the 
results from this routine are shown in Fig. 30 
together with the corresponding data for po{Q^) 
that indicate the ranges selected for Q i and Qt.lt 
is seen that the conditional distributions from this 
simulation show the same stochastic trends seen 
for the experimental data in Fig. 22, i.e., 
(<2 "*■ T=> <^ 7 i) for all values of /. 



The routine shown in Table 6 is used to analyze 
simulated PD pulses in real time. The value for Q * 
is stored from the previous half-cycle. Each time a 
negative pulse is generated, its number, phase, and 
amplitude are tested and selected for inclusion in 
various bins associated with distributions such as 
considered in Table 6. Data for many different 
distributions can be simultaneously recorded in this 
way. Thus it is possible using software to make 
efficient use of available data without generating 
large intermediate data files. 

In order to implement this approach in a 
measurement system, it is necessary that the data 
be converted to digitized pairs of numbers corre- 
sponding to pulse amplitude and time (or phase) 
that can be accessed sequentially by a computer. 
Such a system could, of course, also be used to 
analyze prerecorded data. Although this approach 
appears feasible using existing analog-to-digital 
conversion methods and computer technology, no 
attempts have been made to develop the required 
hardware. 

It should finally be noted that Monte Carlo 
simulations such as the one that yielded the results 
shown in Fig. 30 could be used to test the perfor- 
mance of stochastic analyzers. This can be accom- 
plished by using a digital-to-analog converter to 
produce a sequence of pulses from the output of 
the computer simulation. The conditional and 
unconditional distributions measured for the 
simulated pulses can be compared with the known 
distributions determined from an analysis made 
using computer routines such as given in Table 6. 
The results from the simulation can also be used to 
provide an indication of the quantities of data 
needed to obtain statistically meaningful results. 
The results shown in Fig. 30 correspond to 10* 
cycles of the excitation voltage. It is clear that the 
quantity of data for /)i((^8 Is*) is close to the min- 
imum needed for determination of a statistically 
significant memory effect. 

The possibilities thus appear to exist for 
constructing efficient computer-based stochastic 
analyzers that can essentially duplicate and extend 
the capabilities of the system described here. There 
are no reasons why such a system could not be 
designed to accept either analog or digital data 
directly from measuring devices or indirectly from 
prerecorded files. The correct performance of the 
system can be verified using reference data with 
known stochastic behavior such as generated with a 
Monte Carlo simulation. 
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Fig. 30. Examples of results from a computer analysis of PD pulses generated using a Monte Carlo simulation. Shown are the 
conditional and unconditional distributions determined from a data sorting routine similar to that listed in Table 6. 
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